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HATIOML  ADVISOEY  COMMITTEE  FOE  AERONAUTICS 


TECHNICAL  MEMORANDUM  1277 


PRESENT  STATE  OF  DEVELOPMENT  IN  NONSTEADY  MOTION 
OF  A  LTFTING  SURFACE  * 
By  P.  Cicala 


A  sumnaaiy  is  given  of  the  principal  resiilts  thus  far  obtained 
from  studies  of  the  nonsteady  motion  of  a  lifting  surface  In  an 
inconrpressihle  fluid;  the  methods  followed  hy  YEtrious  investigators 
are  indicated* 

The  aerodynamic  prohlem  of  the  nonsteady  motion  of  an  airfoil 
has  heen  the  suhject  of  numerous  investigations,  vhich  in  20  years 
have  "brought  a  degree  of  development  such  that  an  entire  hranch 
(that  of  two-dimensional  motion)  may  he  said  to  have  heen  completely 
solved.  The  mass  of  existing  puhlications  is  very  large  and  among 
these  naturally  many  overlap;  moreover,  the  study  of  the  prohlem  has 
produced  a  variety  of  methods  such  that  the  same  phenomenon  is 
endowed  with  rather  diverse  aspects  that,  although  of  considerahle 
speculative  interest,  do  not  always  facilitate  the  task  of  those  who 
wish  to  learn  only  the  results  of  the  research.  A  synthetic  deduc- 
tion of  the  results  thus  far  obtained  is  given  herein!  in  order  to 
expound  the  principles  of  the  various  methods  of  investigation  and 
particularly  to  collect  the  latest  results  in  a  form  that  is  heet 
suited  to  application. 

As  previously  stated,  the  two-dimensional  problem  has  heen 
exhaustively  studied.  In  this  field,  the  problems  of  unsteady  motion 
reduce  to  computations  that  at  times  are  laborious,  but  which  can 
always  be  conducted  without  uncertainties  of  the  approximations  and 
which  are  considerably  facilitated  by  existing  tabulations.  Only  a 
few  attempts  have  been  made  to  develop  a  rigorous  computation  for 
the  wing  of  finite  aspect  ratio  inasmuch  as  the  existing  methods 
present  approximate  solutions  that  contain  many  inaccuracies. 

In  part  I  of  this  report,  the  results  relative  to  the  wing  of 
infinite  aspect  ratio  are  described,  always  with  the  assumption  of  a 
perfect  and  incompressible  fluid  with  regard  to  the  components  of 


*"Lo  State  Attuale  delle  Ricerche  sul  Moto  Instazionario  di  Una 
Superficle  Portante."  Estratto  Da  "L'Aerotecnica"  Vol.  XXI,  N.  9-10, 
Settembre-Ottobre  1941,  XIX. 

-'-The  analytic  developments  have  "been  shortened  and  the  rigor  of 
the  demonstrations  has  at  times  "been  relaxed. 
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of  tl^e  normal  force  on  the  wing  plane.  In  peart  11^  procedures  are 
indicated  for  the  computation  of  the  wing  of  finite  aspect  ratio; 
a  compsirlson  is  made  of  the  various  approximations  thus  far  employed. 
In  part  III,  the  results  are  given  relative  to  the  force  component 
in  the  direction  of  motion  (propulsive  or  drag).  In  part  IV,  the 
results  of  experimental  investigations  are  considered.  In  each  part 
em  evaluation  is  given  of  the  existing  publications,  which  may  sei^e 
as  a  guide  to  a  more  detailed  study  of  aiiy  aspect  of  the  prohlem. 

Once  a  reference  system  is  established  with  respect  to  which 
the  fluid  is  motionless  at  infinity,  it  is  assumed  that  a  fixed 
plane  exists  from  which  the  points  of  the  lifting  surface  are  at 
distances  that  may  always  he  considered  very  small  with  respect  to 
the  dimensions  of  the  projection  of  the  surface  on  the  plane.  The 
thickness  and  the  curvature  of  the  profile  are  therefore  considered 
infinitesimal,  as  are  the  displacements  along  the  normal  to  the 
fixed  plane.  Almost  all  the  investigations  up  to  the  present  time 
on  the  nonsteady  motion  of  a  wing  use  simplifications  that  are 
derived  from  the  preceding  assumption  and  thus  assime  that  the  sin- 
gularities in  which  the  "body  in  motion  is  schematized  are  permanently 
contained  in  the  fundamental  plane.  In  this  plane  lie  the  orthogonal 
axes  I  and  y,  which  are  displaced  with  respect  to  the  fixed  refer- 
ence and  remain  parallel  to  themselves  with  a  velocity  V  peirallel 
to  the  X-axis  hut  opposite  in  direction;  the  z-axis  is  at  right  ajigles 
to  the  x,y-axes.  With  respect  to  these  axes,  which  follow  the  motion 
of  the  wing,  the  relative  velocities  of  the  points  of  the  wing  are 
small  with  respect  to  V.  It  is  thus  assumed  that  the  perturbations 
produced  by  the  motion  of  the  wing  are  siifficiently  small.  Hence,  in 
the  relations  that  are  used,  all  terms  of  higher  order  than  the  linear 
terms  in  the  velocities  induced  by  the  motion  of  the  obstacle  in  the 
surrounding  fluid  are  neglected.  This  fundamental  simplification 
gives  the  problem  under  consideration  the  advantages  of  the  linear 
theories  (the  most  important  of  which  is  the  principle  of  superpo- 
sition) from  the  effects  of  which  by  the  analysis  of  particularly 
siniple  motions,  solutions  caji  be  obtained  that  through  linear  combi- 
nations make  possible  the  study  of  motions  of  a  more  complicated 
character. 

The  velocity  V  is  assumed  to  be  small  compared  with  the  velo- 
city of  sound  althoxigh  investigations  have  been  conducted  that  consider 
the  compressibility  of  the  medium  (reference  l). 

The  variations  in  time  or  space  must  frequently  be  measiared. 
The  symbols  occurring  in  the  derivations  herein  are  therefore  defined 
in  the  following  paragraph: 
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Let  q.  be  any  magnitude  that  is  a  fvmction  of  a  point  and  that 
for  every  point  varies  in  time.  The  symbol  ^/dx  (or  ^/^y) 
denotes  the  dependence  of  the  law  of  variation  of  q.  as  a  function 
of  X  (or  of  y),  measured  hy  giving  to  q,  at  each  point,  the 
values  that  correspond  to  a  certain  fixed  time  t.  The  symbol 
dq./dt  indicates  the  rate  of  variation  of  q.  at  a  fixed  point 
relative  to  the  x,  y,  eind  z  axes,  which  in  the  absolute  refer- 
ence system  is  displaced  with  the  velocity  T.  In  general,  with 
the  assimied  linearization,  the  derivative  also  represents  the  rate 
of  variation  of  q.  at  a  fixed  point  of  the  wing  because  the  velocity 
of  a  point  of  the  wing  with  respect  to  the  moving  axes  is  very  smallj 
if  the  gradient  of  q  is  not  too  great,  the  variation  of  q  in  time 
will  be  the  same  whether  measured  at  the  fixed  point  relative  to  the 
moving  axes  or  measured  at  the  point  ttiat  follows  the  lifting  surface. 
In  order  to  express  the  rate  of  change  in  time  of  a  q  always 
measured  at  the  same  point  of  the  stationary  reference,  the  symbol 
d'q^/dt  is  used.  Because  of  the  assumed  linearization,  the  previous 
derivative  coincides  with  the  derivative  that  measures  the  change  in 
q  with  time  for  the  same  material  point.  Inasmuch  as  the  absolute 
velocities  of  the  fluid  particles  are  very  small,  if  the  gradient  of 
q  is  not  too  large,  the  change  of  q  with  time  will  be  the  same 
whether  measured  at  the  fixed  point  or  following  the  molecule.  It  is 
preferred  in  this  report  to  use  a  distinctive  sign  in  d'q  to  differ- 
entiate the  local  derivative  from  the  derivative  of  the  quantities 
that  depend  for  the  problem  under  considei^tion  only  on  the  time  param- 
eter; for  this  problem  the  notation  dq/dt  is  used. 


PAET  I.  -  TWO-DIMEHSIOML  PROBLEM 

1.  The  results  of  modem  research  on  the  two-dimensional  problem 
will  first  be  described.  The  simplification  introduced  by  the  assump- 
tion that  the  phenomenon  develops  in  the  plane  of  the  x-  and  z-axes 
is  such  that  it  can  be  stated  that  each  problem  within  this  range 
can  be  reduced  to  the  computation  of  integrals  that,  with  graphical 
procedures  aided  by  analytical  considerations,  can  be  computed  with 
suitable  accuracy  (limited  to  the  approximation  of  the  existing 
tables),  which  eliminates  the  expsmsion  into  series.  In  the  three- 
dimensional  case  (the  wing  of  finite  span),  except  for  some  parti - 
ular  problems,  only  approximate  solutions  exist. 

In  order  to  simplify  the  expressions  of  the  two-dimensional 
problem,  the  semichord  of  the  wing  is  assumed  to  be  of  unit  length. 
The  abscissas  of  the  leading  and  trailing  edge  of  the  wing  are  assumed 
to  be  given  by  x  =  -1  and  x  «  1,  respectively.  In  order  to  define 
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the  points  of  the  profile,  the  parameter  i^  is  introduced.  The 
variation  with  x  is  shown  "by  the  relation 


cos 


The  values  0  and  k  of  the  parameter  therefore  correspond  to 
the  leading  and  trailing  edges  of  the  wing,  respectively. 

The  coordinate  z,   normal  to  x,   is  considered  positive  in  the 
downward  direction  (fig.  1).  The  vertical  component  w  of  the  velo- 
city of  the  fluid  is  therefore  considered  positive  if  tiimed  in  the 
direction  of  positive  z.  The  difference  in  pressure  p  "between  the 
two  faces  of  the  wing  has  the  positive  sign  directed  upward.  The 
same  siga  convention  is  true  for  the  vertical  force  P,  the  lift  of 
a  segment  of  unit  chord.  The  moment  M  on  the  same  segment  is  con- 
sidered positive  if  it  is  a  diving  moment. 

2.  Condition  of  tangency.  -  The  condition  that  the  profile  he 
impenetrable  to  the  fluid  is  expressed  hy  making  the  relative  velo- 
city of  the  fluid,  with  respect  to  the  wing,  tangent  to  the  profile; 
or,  in  other  words,  the  absolute  velocities  V©  of  a  point  on  the 
contour  and  Vf  of  the  fluid  particle  in  contact  with  the  contour 
have  the  same  projection  on  the  normal  n  (fig*  2).  Inasmuch  as 
only  the  linear  terms  in  the  coordinate  z  of  the  points  of  the  con- 
tour and  their  derivatives,  or  in  the  velocities  (perturbations) 
created  in  the  fluid  by  the  motion  of  the  wing  are  considered,  the 
component  on  n  of  the  velocity  V^  may  be  supposed  eq^ual  to  the 

component  w  of  the  velocity  of  the  fluid  peirticle,  which  is  parallel 
to  z.  By  using  the  simplification  of  neglecting  the  quadratic  terms 
in  the  computation  of  the  projection  on  n  of  the  velocity  V^,, 

which  has  the  components  V  and  ^z/dt,  there  is  obtained 

5.  Circulntion  and  pressure  on  profile  in  steady  motion.  -  Under 
the  conditions  of  steady  motion,  it  is  known  that  the  relative  velo- 
cities of  the  fluid  with  respect  to  the  profile,  even  if  the  profile 
is  considered  to  be  of  infinitesimal  thickness,  axe,  in  general, 
different  on  the  two  surfaces.  In  the  motion  imder  consideration, 
if  Vi  and  V2  are  the  velocities  at  the  corresponding  points  of 
the  two  surfaces  (fig.  1),  the  pressure  rise  p  between  the  surfaces 
is  given  by  the  Bernoulli  equation 
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=  i  p(t2„-v2  )  =  p(V2-Vl)  -^ 


PVu 


•where  it  may  "be  assumed,  "because  of  the  linearization  hypothesis,  that 
the  average  of  the  velocities  Y±    and  V2  is  V,  and  u  denotes 
the  difference  hetween  them. 

If  the  wing  is  in  the  positive  aspect,  the  smallest  velocity  is 
found  at  the  tottom  surface.  The  velocity  difference  u  is  con- 
sidered positive  under  these  conditions  and  the  positive  pressure  p 
is  therefore  in  the  opposite  direction  to  positive  z.  la  order  to 
represent  the  field  created  hy  the  wing,  the  skeleton  of  the  wing 
composed  of  a  vortex  film  is  considered;  this  system  of  singularities 
is  capahle  of  giving  the  existing  velocity  increment  u  "between  the 
two  surfaces  if  the  circiilation  in  an  element  dx  is  equal  to  u  dx. 

4.  Circulation  and  pressure  on  profile  in  unsteady  motion.  - 
For  the  case  where  the  motion  is  unsteady,  the  velocity  increment  may 
"be  represented  by  a  vortex  distri"bution  of  Intensity  g  «  u.  The  only 
difference,  when  compared  with  the  preceding  case,  is  in  the  fact  that 
the  discontinuity  in  the  velocity  field  exists  not  only  on  the  points 
of  the  wing  "but  also  in  the  wake  "behind  the  wing.  In  order  to  descri"be 
this  phenomenon,  consider  two  fluid  layers  that  pass  above  6uid  below  the 
profile.  Because  of  the  linearization  assumption,  the  distances  of  the 
points  of  the  profile  from  the  x-axis  can  be  neglected  in  the  following 
discussion.  It  can  easily  be  verified  that  the  terms  of  the  second 
order  will  therefore  be  neglected.  The  difference  between  the  momen- 
tums  of  the  two  layers  is  computed  (fig.  3)^  at  the  time  tg  and  at 

the  time  tg  +  dt  and  therefore  the  variation  that  the  difference  has 

undergone  in  the  interval  dt  is  also  computed.  This  variation, 
divided  by  the  interval  dt,  mnst  be  equal  to  the  difference  between 
the  forces  along  x  that  are  applied  to  the  two  layers.  Inasmuch  as 
tangential  actions  do  not  exist  (perfect  fluid),  and  becaxise  no  pressure 
difference  exists  on  the  anterior  face  (upstream  of  the  profile),  the 
previously  mentioned  momentums  will  be  given  by  (pi-P2)  ^Ly  «  p  dy. 
The  difference  in  the  momentums  of  the  molecules  of  the  two  layers, 
which  is  given  by  the  product  of  the  mass  and  the  difference  in  the 
velocities  u,  varies  in  the  interval  dt  because: 

2 
The  two  layers  are  Indicated  in  figure  3  by  hatched  lines  in  the 

two  directions;  the  position  of  the  leading  edge  of  the  profile  is 

indicated  by  a  semicircle. 
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1.  New  molecviles  come  in  contact  with  the  profile^  and  there- 
fore aGq.mre  the  velocity  increment  u  =  g. 

2.  The  increment  for  a  given  point  of  the  profile  varies  with 
time. 

If  the  graphs  of  g  relative  to  t  =  to  and  t  =  tg  +  dt  on 

the  same  position  of  the  profile  (fig.  3)  are  plotted  corresponding 
to  the  two  causes,  respectively,  it  is  found  that  the  area  (which 
gives  the  difference  in  the  momentums  when  multiplied  hy  P  dy)  varies 
in  the  interval  dt  hy  the  amount 


gV  dt  +  dt 


a, 
-1 


t^- 


The  first  term  is  indicated  in  figure  3  by  the  obliquely  hatched 
area  (which,  except  for  infinitesimals  of  higher  order,  is  equivalent 
to  a  rectangle  of  base  V  dt  and  altitude  equal  to  the  value  of  g 
corresponding  to  the  abscissa  x  of  point  P);  the  second  term  is 
indicated  by  the  vertically  hatched  area.  Equating  the  impulse  to 
the  variation  of  the  momentum  and  dividing  by  p  dy  dt  yleld^ 


^     -J 


-1 


^  dx  (2) 


The  first  term  on  the  right  side  measures  the  pressure  at  P, 
which  is  obtained  under  the  conditions  of  steady  motion;  this  pres- 
sure depends  on  the  local  value  of  the  velocity  increment.  By  the 
effect  of  tJie  other  term,  the  pressure  under  the  conditions  of  non- 
steady  motion  depends  on  the  variation  that  g  undergoes  at  the 
instant  considered  In  the  entire  strip  ahead  of  P. 


"The  difference  of  the  momentums  of  the  two  layers  is  observed  to 
be  zero  befors  arriving  at  the  wing  because  the  velocity  Increment  in 
front  of  the  wing  is  zero  if  there  are  no  other  lifting  surfaces. 

This  relation  is  obtained  by  another  method  in  differential 
form  in  reference  2, 
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5.  Yortlclty  of  wake-  -  If  at  the  instant  t  =  to  the  point  P 
coincides  with  the  trailing  edge  of  the  wing,  equation  (2)  still 
holds,  "but  the  first  term  becomes  zero  hecause  a  difference  in  pres- 
sure cajmot  exist  where  there  is  no  wing  surface.  It  follows  that 
the  increase  in  velocity  g„  existing  in  the  points  of  the  wake, 

which  is  zero  in  the  case  in  which  the  distribution  of  g  on  the 
wing  is  independent  of  the  time,  is  given  by 

^  (3) 


V  dt 


where 


■i: 


K  =     g  dx 


is  the  circulation  about  the  wing  and  the  derivative  is  measured  at 
the  instant  at  which  the  trailing  edge  passes  through  the  point  of 
the  wake  under  consideration.  In  other  words,  in  the  distance  V  dt 
that  the  wing  moved  in  the  time  interval  dt  a  vorticity  is  distri- 
buted equal  to  the  variation  that  the  circulation  about  the  profile 
has  undergone  in  the  same  time.  This  conclusion  can  also  be  derived 
from  the  principle  of  the  conservation  of  vorticity. 

In  figure  4,  the  diagram  of  the  vorticitles  on  the  profile  axid 
in  the  wake  is  given  for  the  case  in  which  the  wing  executes  a  trans - 
latory  oscillation  with  frequency  Q   related  to  L  and  V  by 

QL  =  2V 


The  curve  a  refers  to  the  instant  in  which  the  wing  crosses 
the  middle  position;  b  refers  to  the  position  at  the  end.  The 
scale  of  the  vortices  is  indicated  by  assuming  the  vertical  maximum 
velocity  of  the  profile  to  be  unity.  For  compajrlson,  figure  4  shows 
the  graph  of  the  circulation  g'   corresponding  to  the  maximum  veloc- 
ity in  the  case  of  steady  motion.  The  increment  of  velocity  in  the 
wake  is  considerable.  In  the  theories  of  a  wing  in  nonsteady  motion, 
it  is  generally  assumed  that  this  velocity  increment  remains  in  the 
position  in  which  it  originated,  or,  in  other  words,  that  the  vortices 
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shed  "by  the  ving  maintain  their  original  position  and  intensity- 
unaltered  in  time.   It  is  evident  that  the  error  inherent  in  this 
assumption,  from  which  the  real  phenomenon  certainly  deviates, 
impairs  to  some  extent  the  results  of  the  theory. 

6.  Bound  and  free  vortices.  -  An  original  physical  interpre- 
tation  of  equation  (2)  was  given  hy  Bimlaaum  (reference  3).  By 
setting 


--I 


3f^  (^) 


Equation  (2)  is  written  as 

g  =  7+  €  (6) 


With  the  aid  of  this  relation,  the  value  of  the  total  vorticity 
g  at  a  point  on  the  profile  at  which  the  velocity  increment  between 
the  corresponding  points  of  the  top  and  bottom  surfaces  is  divided 
into  two  parts : 

(1)  The  bound  vorticity  7,  which  is  sustained  by  the  aerody- 

namic action 

(2)  The  free  vorticity  £,  which  trails  in  the  fluid  in  its 

relative  motion  and  which  therefore  gives  no  pressure 
rise 

In  order  to  clarify  the  relation  between  the  free  vortices  and  the 
bound  vortices,  equation  (2)  is  differentiated  with  respect  to  x 
after  g  is  expressed  in  terms  of  eqviation  (6).  Thus, 


¥.- 


3£   „  de   d'€  ,^> 

5t^^57  =  -dt  (') 
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This  relation  indicates  that  the  variation  undergone  in  an  inter- 
val of  time  by  the  vortex   e  at  a  point  of  the  fluid  is  equal  (and  of 
opposite  sign)  to  the  variation  undergone  in  the  same  time  "by  the 
bound  vortex  at  the  point  of  the  profile  in  contact  with  it.  In 
other  words,  the  bound  vortices  leave  at  each  instant,  in  the  fluid 
with  which  they  come  in  contact,  an  effect  represented  by  a  vorticity 
of  intensity  equal  to  that  which  they  have  produced  at  that  instant. 
According  to  this  representation,  each  of  the  inducing  elements  at 
the  wing  is  considered  in  isolation  and  therefore  produces  a  pres- 
sure rise  expressed  by  equation  (4)  and  has  a  proper  vortex  wake. 
The  free  circulation  at  any  point  on  the  wing  or  downstream  of  it 
is  given  by  the  sum  of  the  circulations  of  the  wakes  corresponding 
to  the  boiHid  vortices  that  are  upstream  of  the  point  considered. 

The  decomposition  of  the  total  circulation  in  the  two  vortex 
systems  previously  described  is  shown  in  figure  4.  It  can  be  seen 
that  the  graph  of  the  free  circulations  on  the  chord  are  joined  in  a 
continuous  manner  with  the  system  of  the  wake  at  the  point  where  the 
free  circulation  is  equal  to  the  total  circulation. 

7.  Relations  between  circulations  and  normal  velocities.  - 
Ordinarily,  the  law  of  variation  of  w  along  the  chord  and  with  time 
is  known  from  equation  (1),  with  the  aid  of  which  these  quantities 
are  derived  from  the  characteristics  of  the  motion.  The  vorticities 
on  the  wing  and  in  the  wake  must  induce  at  each  Instant  on  the  points 
of  the  chord  the  assigned  w;  that  is. 


=   2«j- 


-M=^l^^^^  (8) 


The  integration  is  extended  from  the  leading  edge  of  the  wing 
to  the  entire  part  of  the  wake  in  which  the  passage  of  the  wing  has 
created  the  vorticity  discontinuity. 

It  is  known  that,  for  the  wing  in  steady  motion,  the  field  of 
motion  would  not  be  determined  if  the  point  of  separation  (trailing 
edge  of  the  wing)  were  not  fixed.  This  condition  is  also  assumed 
for  the  wing  in  unsteady  motion  and  is  translated  into  the  analytical 
condition  that  g  be  finite  at  the  trailing  edge  of  the  wing. 

The  vorticity  in  the  wake  is  connected  with  that  at  the  wing 
by  equation  (3).  When  equation  (8),  which  is  completed  by  the  condi- 
tion of  separation  and  by  equation  (3),  is  solved  the  total  vorticity 
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g  is  oTDtained.  The  pressures  are  then  obtained  "by  making  use  of 
equation  (2).  As  vrill  he  shown  in  more  detail  in  section  17^  this 
scheme  is  followed  by  various  procedures  of  solution  hy  Theodorsen, 
who  makes  use  of  conformal  transformations  (reference  4),  Schwarz 
(reference  5),  who  makes  use  of  Betz's  solution  of  equation  (8),  and 
Sohngen  (reference  6).  Wagner  and  Glauert  also  refer  to  the  total 
vorticity,  but  determine  the  total  action  on  the  profile  and  not  the 
pressure  distribution  (as  do  von  Ka'rma'n  and  Sears,  reference  7). 

When  a  different  method  is  used,  the  free  circulation  on  the 
chord  and  in  the  wake  can  be  expressed  as  a  function  of  7  by 
nakiag  use  of  the  integrated  equation  (7)  and  thus  transforming 
equation  (8)  so  that  only  the  unknown  function  7  appears  in  it. 
Then  the  pressures  can  be  directly  obtained  by  means  of  equation  (4). 
This  procedure  was  followed  by  Birnbaum  (reference  3),  by  Kussner 
(reference  8),  and  by  Cicala  (reference  9). 

8; Acceleration  potential.  -  A  different  interpretation  of  the 
same  problem  can  be  made  on  the  basis  of  the  acceleration  potential. 
For  a  perfect  and  incompressible  fluid,  the  equation  of  Euler 


a  =  -  grad  £ 

expresses  the  equality  between  the  force  of  inertia,  which  corresponds 
to  the  acceleration  a  of  the  fluid  particle,  and  the  resiiltant  of 
the  forces  that  are  transmitted  to  the  particle  by  the  medium  surround- 
ing the  particle.   The  Euler  equation  also  permits  stating  that  the 
components  of  the  acceleration  can  be  obtained  from  the  fimction  p/p 
through  the  same  operations  of  differentiation  with  which  the  com- 
ponents of  the  velocity  are  deduced  from  the  corresponding  potential. 
The  generating  singularities  are  arranged  on  the  surface  of  discontinu- 
ity (wing  +  wake)  for  the  acceleration  field  whose  potential  satisfies 
analytical  properties  similar  to  those  of  the  velocity  potential  in 
the  same  manner  as  for  the  velocity  field.  The  singularities  are 
arranged  where  the  discontinuity  exists  in  the  pressures;  that  is, 
on  the  lifting  surface.  The  simplification  that  is  introduced  by 
this  concept  is  not,  however,  as  great  as  might  appear  from  the  fact 
that  the  singularities  (which  are  called  the  dipoles  of  the  pres- 
s-'re)  are  limited  only  to  the  wing.  In  each  case,  it  is  necessary 
to  pass  from  the  acceleration  field  to  the  velocity  field  and  this 
pa=:sage  requires  an  integration  through  which  the  effects  of  all  the 
preceding  states  of  motion  are  felt  and  which,  according  to  the 
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vortical  representation,  leaves  a  trace  in  the  system  of  free  vortices- 
In  other  words,  the  velocity  field  of  the  pressure  dipoles  depends 
on  the  history  of  the  formation  of  the  dlpole. 

The  relation  TDetween  the  vortex  representation  of  the  phenomenon 
and  that  "based  on  the  pressure  dipoles  is  already  implied  in  the  con- 
cept of  BirnlDaum  of  the  vortices  shed  from  all  the  points  of  the 
airfoil  that  is  the  only  source  of  the  hound  vorticity  7,  which  is 
proportional  to  the  difference  in  pressiore  between  the  two  faces  of 
the  profile.  A  method  will  be  shown,  on  the  basis  of  the  concept  of 
vortices,  for  the  derivation  of  some  fundamental  relations  that  in 
other  publications  are  justified  with  the  procedure  of  the  accelera- 
tion potential. 

The  system  of  bound  vortices  induces  at  a  point  0  at  the 
instant  t  =  ^  a  velocity  that  is  denoted  Vg^;  at  a  succeeding 

instant  t  =  t  +  a ''",  at  a  point  that  occupies,  with  respect  to  the 
wing,  the  same  position  that  0  occupied  in  the  first  condition, 
there  will  be  a  velocity  differing  from  Vg^  by  an  amount  denoted 

by  AVg^.  The  velocity  AV^,   will  depend  on  the  variation  in  the  con- 
stitution of  the  induction  system  or,  with  changed  sign,  will  be  the 
velocity  Induced  by  the  system  of  vortices  that,  at  the  instant  con- 
sidered, are  freed  from  the  bound  system  and  constitute  the  trace 
that  the  boujid  system  has  left  in  the  fluid  with  which  it  has  been 

in  contact.  The  quantity  AV  /a T  ,  as  A T-yQ  as  a  limit,  represents 

a 

the  derivative  ^V^/dt-     Hence,  if  it  is  desired  to  express  analyti- 
cally the  property  that  the  actual  velocity  Yf  of  the  fluid  parti- 
cles results  from  the  sum  of  the  velocities  induced  by  the  actual 
configuration  of  the  bound  vortices  and  from  those  freed  at  the  pre- 
ceding instants,  then 


^t 


Yf  =  Ya 


where  Va  a-^-d  Vf  are  computed  at  tne  point  and  at  the  instant  t, 
3v/(3t  being  computed  at  the  point  of  the  fixed  space  being  considered, 
and  the  system  of  boimd  vortices  is  in  the  corresponding  position  at 
the  instant   t  preceding  t.   In  the  case  of  translational  motion 
or,  more  precisely,  in  the  linearization  assumed,  the  terms  of  the 
second  order  in  the  velocities  of  the  points  of  the  wing  relative  to 
a  system  that  is  displaced  with  translational  motion  with  velocity  V 
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are  neglected;  the  term  dVa/^t  is  measured  by  having  the  point 
of  Induction  at  the  ahBCissa  x-V(t-T),  or,  with  respect  to  the 
point  0  of  the  ahscissa  x  to  which  the  velocity  Vj  refers, 
in  a  position  upstream,  advanced  hy  the  quantity  V(t-T)   if  V  is 
assumed  constant. 

The  following  relation  can  he  obtained  from  equation  (9): 

ilS:  =  ^  V  (10) 

dt    SjT 

The  quantity  Vg^  is  the  velocity  induced  at  0  at  the  constant 
t  by  the  bound -vortex  film  layer.   From  the  instant  t-dt  to  the 
instant  t,  the  vortex  system  changes  only  in  the  number  of  the  bound 
vortices  leaving  a  trace  and  that  occupy  the  position  corresponding 
to  the  time  t.  The  variation  that  the  velocity  Vf  of  the  fluid 
has  undergone  at  the  point  0  of  the  fixed  space  because  of  the 
effect  of  this  change  (local  variation  d'Vj>/dt  dt)  is  that  which 

would  be  measured  by  leaving  the  vortices  7  in  position  and  moving 
the  point  of  induction  of  the  segment  dx  =  Vdt  so  that  the  point 
of  induction  occupies  the  position  that  corresponds  to  the  instant  t; 
the  variation  that  is  measured  in  this  manner  is  expressed  by 
dVa/^x  dx  and  is  thus  added  to  eqioation  (10),  which  equates  the  two 
expressions  of  the  variation. 

Relations  (9)  and  (10)  are  equivalent  to  the  relations  derived 
by  Possio,  which  are  based  on  the  concept  of  the  acceleration  poten- 
tial, if  it  is  assumed  that  the  stationary  field  C,  defined  in 
referenc*te  10,  coincides  with  the  field  produced  by  the  system  of 
bound  vortices  in  a  uniform  stream  V-   Inasmuch  as  the  concept  of 

Equation  (lO)  is  derived  directly  If  the  acceleration  field  is 
considered  to  depend  only  on  the  actual  values  of  the  pressures  on  the 
profile  and,  with  the  pressures  equal,  has  the  same  configuration  as 
though  the  motion  were  steady-   In  this  case,  according  to  equa- 
tion (4),  the  actual  pressiires  would  be  obtained  by  having  only  the 
actual  values  of  7   on  the  profile:   the  velocity  of  the  fluid  would 
be  Va  and  the  acceleration  that,  by  the  assumed  linearization,  is 
computed  from  d'Va/dt  (and  not  from  the  derivative  formed  in  following 
the  fluid  particle)  would  be  expressed  by  d'Va/dt  =  djei/dx  V   (without 
dVa/<3t  because  this  virtual  field  is  steady).  When  the  acceleration 
of  this  virtual  field  is  equated  to  the  effective  acceleration,  which 
is  written  d'Vf/dt,  equation  (10)  is  obtained,  which  when  integrated, 
gives  equation  (9). 
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the  acceleration  potential  lends  itself  less  to  physical  interpre- 
tations than  the  concept  founded  on  vortices  although  some  simpli- 
fications can  he  obtained  in  the  analytical  development^  it  is 
preferred,  in  the  following  sections,  to  employ  the  classical  method 
of  description. 

9.  Cases  of  total  zero  circulation-  -  If  the  total  circulation 


ahout  the  wing  is  constant  in  time  and  the  system  has  no  vortices  in 
the  wake,  the  induced  velocities  can  be  computed  on  the  basis  of  the 
total  circulation  about  the  wing  at  the  instant  considered  in  the 
same  manner  as  for  steady  motion.  From  the  analysis  relative  to  this 
case  (reference  11,  p.  185),  it  is  known  that  if  the  circulations  are 
represented  by  one  of  the  following  functions 


Si 


cot 


^ 


2  sin  i» 


g^^  =  -  2  sin  n^^  (n  =  2,  3,  ...) 


(11) 


the  corresponding  velocities  are  respectively  given  by 


wj_  =  2  "^  cos  i> 


v^   =  cos  n^ 


(12) 


Inasmuch  as  the  values  of  g^  represented  by  equation  (11) 
itisfy  the  condition 


K 


dx  =  0 


the  preceding  result  holds  for  any  motion  because  no  wake  exists  down- 
stream of  the  body.  Thus, 


w  =ZAn  Wn 
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where     An     are  qioantities  independent  of    iJ      but  fionctions  of   time. 
The  total  circulations  are  represented  by  the  corresponding  sum 


g  =  ^An  gn 


(13) 


The  corresponding  hound  vortices  are  represented  hy  S  y^,   where, 
as  derived  from  equations  (5)  and  (6), 


'n 


7-,  =  A-j  (  cot  p  -  2  sin  ^  )  +  A'-^  (sin  <>  +   sin  i>  cos  ^) 


>(14) 


^^   V  dt 


(15) 


In  general,  the  values  of  w  can  he  expressed  hy  expanding  in 
a  Fourier  series: 


W  =  —  Ag  +  2Aj^  cos  Xi'i^ 


(16) 


The  coefficients  An,  general  functions  of  time,  can  he  obtained  by 
harmonic  analysis  by  setting 

w  =  w  -  lAn  wn 

There  is  easily  obtained 


w  =  2  (Ao 


-  Ai) 


(17) 


It  is  therefore  concluded  that,  if  the  vertical  velocities  on 
the  profile  are  developed  in  Fourier  series,  the  circulations  and 
the  pressures  can  be  directly  computed  with  the  aid  of  equations  (13) 
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and  (14)  if  Aq  =  A^.   In  the  general  case,  It  will  "be  necessary  to 
sum  the  circulations  (eqioation  (13)),  which  correspond  to  the 
velocity 


-  1  r" 

w  =  —  1   w 


(1  -  cos  ■^)  di>  (17'  ) 


which  is  a  function  of  time,  but  is  constant  over  the  chord. 

10.   Case  of  velocity  v  constant  over  chord-  -  Applying  equa- 
tion (10)  for  the  potential  on  the  z-axis,  inasmuch  as  w  is  the 
projection  of  V.^,  yields 

d'v    ^a  (      . 

vlt  =  Sr  ^^®^ 

where  Va  is  the  projection  of  the  vector  Y^     on  z. 

Because  at  each  instant  w  assumes,  for  all  the  points  of  the 
chord,  the  value  w,  the  following  relation  is  obtained : 


d'w  _  dw 
dt   dt 

and,  therefore,  from  equation  (18),  integrating  along  the  chord  yields 

Tr     dw      _ 
^a  =  ^Tdt  ^  ^  ^ 

where  C  is  a  constant  with  respect  to  x,  hut  is  a  fimction  of  time. 
From  the  analysis  of  steady  motion,  it  is  known  that  the  distribution 
of  the  vorticity  capable  of  giving  velocities  satisfying  the  condition 
at  the  edge  is  given  by 


^0   =  2  ^-^   siniS  +  2  0  cot  ^  (19) 
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In  order  to  investigate  the  dependence  of  C  on  the  condition 
of  motion,  the  elementary  case  is  considered  in  which  v,   ^ich  has 
been  zero  for  an  indeterminate  time,  suddenly  acquires,  at  the 
instant  t',  the  value  Aw  and  maintains  it  unchanged.  For  such 
a  case,  for  t  >  t' ,  C,  which  must  he  proportional  to  Aw,  may  "be 
put  in  the  form 

C  =  (1-E)  Aw  (20) 

where  E  is  a  function  of  the  space  passed  through  by  the  wing 
since  the  instant  t' .  Because,  with  the  passing  of  time,  the  phe- 
nomenon tends  toward  the  steady  conditions  for  which  the  circula- 
tion tends  to  assume  the  distribution 

7q  =  2  Aw  cot  -  (21) 

the  asymptotic  value  of  R  must  be  zero.  The  law  of  the  variation 
of  E  was  studied  in  one  of  the  first  publications  on  the  wing  in 
unsteady  motion  (reference  12,  which  gives  a  resume  of  the  work  of 
Pistolesi  to  whom  reference  is  made).  The  more  general  case  is 
obtained  from  the  elementary  case  by  superposition  of  the  effects. 
According  to  equation  (20),  the  second  term  of  expression  (19)  is 
decomposed  into  the  asymptotic  term  (equation  (21))  and  a  term  that 
contains  the  function  R  and  represents  the  distribution  of  the 
circulation  that  would  be  realized  if  w  possessed,  from  an  indeter- 
minate time,  the  value  Aw  up  to  the  instant  t',  and  then  for 
t  >  t' ,  w  =  0.  The  pressures  corresponding  to  this  term,  which 
represents  the  effect  of  the  preceding  variation  and  diminishes  to 
zero  with  time,  can  be  referred  to  as  "transitory  pressure." 

In  the  computation  of  C  in  a  general  case,  summing  the  effects 
of  all  the  increments  that  w  has  received  from  the  £tart  of  the 
motion,  when  it  may  be  supposed  w  =  0  (so  that   I  Aw  =  w),  there  is 
obtained 


C  =  w 


"w 

E  dw  (22) 

0 
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The  integral  is  taken  as  the  Bum  of  the  products  of  all  the  vari- 
ations that  V  has  undergone  in  the  preceding  instants  (the  procedure 
is  not  modified  -whether  the  variations  are  abrupt  or  gradual)  for  the 
corresponding  values  of  R.  The  distance  traveled  hy  the  wing, 
measured  from  an  arbitrary  origin,  is  denoted  "by  sj  the  value  of  s 
for  the  position  at  which  the  pressures  are  measured  is  denoted  by 
Sq  and  the  value  of  s  for  the  distajace  Sq-s  referred  to  the  seml- 

chord  is  denoted  \iy    a-      The  values  of  the  function  R(o)   are  given  in 
table  I.  These  values  have  been  obtained  from  the  recent  tabulations 
of  Kussner  and  Schwarz  (reference  13). 

11.  General  solution.  -  When  the  ezpresslons  of  the  two  pre- 
ceding sections  are  collected,  the  resultant  pressures  for  a  general 
case  of  motion  may  be  computed.  If  the  values  of  w  are  expressed 
by  the  series  of  equation  (16),  which  is  put  in  the  form 


w  s"  w  +  SAq  Wn 


(16a) 


the  pressures  corresponding  to  the  first  term,  according  to  eqva.- 
tions  (19)  and  (22),  are  given  by 


p  =  PV 


dw 


^rit  ^^^-^H"^ 


••-') 


E  dw  /  cot  ;r 
0        /      2 


whereas,  the  values  of  w  that  constitute  the  summation  of  eqya.- 
tion  (16a)  correspond  to  the  pressures 


p  - PT  Lr 


n 


vhere     y^     is  expressed  by  equation  (14). 

It  is  convenient  to  divide  the  total  pressures  into  a  part  that 
depends  on  the  history  of  the  motion  and  has  been  denoted  as  the 
transitory  pressure  given  by 


PV 


2  cot  I 


E  dw 


(23) 
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and  into  a  part  f,  vhlch  depends  only  on  the  actual  values  of  the 
parameters  that  cheiracterlze  the  motion  and  is  denoted  hy  the 
"instantaneous"  pressure,  given  hy 

P     _     ■^     dw     ^ 
py  =  2  w  cot  2  +  2  y-^  sin  i>  +  Z 7n 

When  equations  (14)  and  (17)  are  used,  this  relation  assumes 
the  form 


^-*o-*|-^  " 


t(^.  *  ''"V"'"-')  ='- » '        '2*) 


For  the  computation  of  the  pressures  on  the  profile,  the  fol- 
lowing operations  are  therefore  required : 

(l)  P^om  the  law  of  motion  there  is  obtained,  with  the  aid 
of  equation  (1),  the  expression  for  w  as  a  function 
of  time  and  of  the  coordiiiate  ^. 

(ll)  By  developing  w  in  a  Fourier  series,  the  values  of  A 
as  functions  of  time  are  computed. 

(Ill)  From  the  instantaneous  values  of  A  and  the  derivatives, 
the  values  of  f     are  computed  from  equation  (24). 
When  w  has  heen  determined  from  equation  (17),  the 
values  of  p  are  computed  hy  equation  (23)  and  hence 
the  resultant  pressure  p  =  p  +  g.  The  total  force 
and  moment  are  obtained  by  the  simple  integrals 


M 


px  dx 


O 


-1 
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The  formulas  given  here  are  also  valid  without  change  if  V 
varies  from  one  instant  to  the  next.  In  such  a  case,  it  is  con- 
venient to  assiome  as  a  reference  variable  instead  of  time  the  dis- 
tance s  traveled  hy  the  wing.  Therefore 

dA^ 

A'n  =  :; — :  where  ds  =  V  dt. 
"   ds 

In  the  calculation  of  the  transitory  pressures,  the  expres- 
sion (23)  is  usually  computed  "by  graphic  integration.  When  the  dia- 
gram of  w  as  a  function  of  0  is  known  and  the  value  of  E  is 
obtained  from  table  I,  it  is  convenient  to  draw  the  graph  of  w(B) 
(fig.  5),  which  is  obtained  by  laying  off,  for  each  ordinate  w, 
the  corresponding  value  of  R  for  the  same  position.  The  area 
enclosed  by  this  curve,  by  the  vertical  axis,  and  by  the  two  hori- 
zontal lines  through  the  ends  (one  of  which  is  the  axis  w  =  0) 
represents  the  integral  of  eq^uatlon  (23)  (cross-hatched  area  in 
fig.  5).  If  when  tracing  around  the  contour  from  a  to  b  the 
area  is  on  the  right,  the  6Lrea  is  considered  positive;  in  figu'^e  5  _ 
the  value  would  be  considered  negative.  Even  though  the  graph  of  w 
presents  abrupt  variations  (as  in  the  case  of  the  figure),  no  compli- 
cations are  thereby  introduced. 

When  it  is  possible  to  proceed  by  the  analytic  method  in  com- 
puting the  preceding  integral,  it  is  convenient  to  assume  as  the 
variable  of  integration  the  distance  s.  Equation  (23)  is  there- 
fore written  as 


£  =  -  3  cot  I 


OSq 


gR(so-s)  ds  (23') 


or  as 


•I 


i.Scot|r°"-gB(0)4  0  (33") 


0 

12.  Example  of  application.  -  Let  it  be  assumed  that  the  wing 
always  displaced  with  constant  velocity  in  magnitude  and  direction 
undergoes  a  sudden  small  rotation  (of  small  amplitude)  about  a  point 
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of  the  chord.  Let  there  "be  determined  the  law  of  the  variation  of 
lift  and  of  the  focal  moment  at  the  instant  of  the  rotation. 

It  is  convenient  first  to  ass\ame  that  the  rotation  occurs  in 
a  finite  interval  of  time  and  then  passes  to  the  limit  to  let  the 
interval  approach  zero.  Let  a  be  the  angle  of  rotation,  cp  the 
final  value  of  a,  and  ^  =  ^q  the  coordinate  of  the  axis  of  rotation. 

While  the  wing  travels  through  the  distance  from  s=0  to  s=A, 
the  angle  a  increases  continuously.  Let 


for      s  <  0,      w  «  0 


for  0<s<A,     w  =  aV  +  T—  (cos  ^o   -  cos    ^) 


and  for   s  >  A,     w  =  ^V 


Hence, 


V^a  +  gcos^o) 


^   „  Z'  _   da 
2 


-1--S 


>?.v(a.r||) 


r  »  ^  +  cos  i^Q 


for  s  between  zero  and  A  the  rotation  is  w  ■  V<p.  In  the  first 
phase. 


ds   \ds     ^3 


f)^ 
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and  in  tiie  second  phase,  dw/ds  =  0.  If  s  =  sq  ("^1^®  position  for 
whicii  the  pressures  are  measured),  let  Kq,   R'q  ...  "be   the  values 
of  E  and  the  derivatives  for  a  =  sg*  Equation  (23' )  is  then  written 
by  expanding  R  in  a  power  series  in  s : 


-p  =  Epv''  cot  I  = 


PA 


(Ro  -  B  R'o  + 


V  /  da     d^a  \  , 
...)  \  -r—  +  r   — 5-  Ids 
^ds     ds2  J 


(25) 


Then 


Jo 

da  , 
ds*^^ 

^A 

^2 

da, 
— 5"  ds 

ds^ 

Jo 

^A 

s 

0 

ds2^^ 

cp 


=  0 


cp 


and  moreover,  the  quantities 


n  da  , 

s     T—  ds 

ds 


s^+1  ^  ds 

ds'^ 


for     A     approaching  zero  "become  zero  at     A^. 
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Thus,  from  ecLuation  (25),  for     A     approaching  zero, 

-p  >  gpcpv"^  cot  _  (Rq   -  r  R'o) 

From  this  relation  it  is  concluded  that  the  transitory  pressure 
decreases  according  to  the  function  R  if  r  «  0  or  if  the  rota- 
tion occurs  ahout  the  neutral  rear  point.  If  r  is  different  from, 
zero,  a  term  is  added  in  the  law  of  vaxiation  of  p  that  decreases 
as  the  derivative  of  R;  this  term  corresponds  to  the  pressure  dis- 
trihution  that  is  created  on  the  wing  in  uniform  rectilinear  motion 
and  that  executes  an  instantaneous  displacement  in  the  direction 
normal  to  the  trajectory  and  then  continues  with  the  initial  speed 
and  direction. 

The  pressures  f  after  the  rotation  axe  given  hy 


f  =  2pcpY^  cot  I 


Hence,  the  lift  after  the  rotation  is  expressed  hy 


I 


(p  +  p)  dx  =  jtpcpLV^  (1  -  R  +  rR*  )  (25d) 


where  L  is  twice  the  chord  of  the  profile  and  R  euid  R'  are  the 
function  of  table  I  and  its  derivative,  respectively,  hoth  approaching 
zero  with  an  increase  of  the  independent  variahle,  which  is  repre- 
sented hy  the  distance  of  the  actual  position  from  that  at  which  the 
rotation  has  occiirred. 

In  take-off  of  the  wing,  that  is,  when  the  wing  starts  its 
motion  from  rest. 

In  the  interval  0  <  s  <  A,  V  passes  from  zero  to  its  final 
value  V„.  For  A— >0, 
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11m 


B  ds 


dv 

ds 


=  Rq  *^C 


and  the  lift  is  therefore  still  expressed  "by  equation  (25a),  in  lirfiich 
the  factor  in  parentheses  reduces  to  1-E  (Wagner's  case).  For  the 
moments,  "because  the  pressures  are  always  proportional  to  cot  i^/2. 


Mp- 


(P  +  P) 


G4) 


dx 


In  the  preceding  computation,  the  impulsive  pressures  that  are 
generated  at  the  instant  of  the  rotation  and  cease  when  rotation 
has  occurred  were  considered.  The  values  are  therefore  immediately 
o"btained  with  the  aid  of  equation  (24). 

15.  Computation  in  finite  terms  of  instantaneous  pressures.  - 
For  the  detennination  of  the  instantaneous  presstires,  a  convenient 
expression  is  given  ty  Sbhngen,  "by  means  of  which  these  expressions 
are  o"btained  directly  from  w  without  expanding  in  a  Fourier  series. 
The  expression,  modified  to  conform  with  the  notation  used  herein, 
may  "be  given  in  the  following  form: 


p/p  Y  «  Aq  cot  ^ 


n« 


-  Bin  ^ 


H(^'.s)  d^' 
cos  ■^'    -  cos  ^ 


(26) 


where 


H(^',s)  -  w(x»,s)  + 


r 


Sw 
3s^ 


dx 


(27) 


find  where  x 
referred  and 


=  -  cos  1^  is  the  coordinate  of  the  point  to  which 
i'  ■  -  cos  ^*     is  the  varla"ble  of  integration. 


p  is 
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The  lower  limit  of  the  integral  in  equation  (27)  is  arhitrary, 
80  that  H  is  defined  except  for  an  additive  constant.  This  arhi- 
trfiLTiness  does  not  affect  the  results  "because 


Off 


d^ 

cos  ^     -  COS  ^, 


0 


0 


0 


The  proof  of  equation  (26)  can  he  given  in  a  manner  that  is  not, 
however,  entirely  satisfactory  from  the  mathematical  point  of  view 
hy  substituting  the  expression  of  equation  (16)  for  w  In  equa- 
tions (26)  and  (27)  and  verifying  that  the  relation  thus  obtained 
agrees  with  equation  (24). 

The  practical  computation  of  the  Integral  that  occurs  in  equa- 
tion (26)  presents  difficulties  for  the  singularity  of  the  function 
integrated  at  the  point  i^  =  ^q.  For  cases  that  are  encountered  in 

practice,  however,  by  dividing  the  chord  into  a  certain  number  of 
strips  it  is  found  that  in  each  strip  H  can  be  represented  by  a 
combination  of  a  few  terms  of  the  type  cos  n^.  Hence,  to  simplify 
the  applications,  Sbhngen  gives  the  following  formula,  which  in  this 
case  permits  conducting  the  computation  in  closed  form: 


o« 


sin  ^ 


f  d^' 


cos  ni 


cos  1^  -  cos  d 


In 


[i- 


cos 


c^- 


cos 


(^-1^2)] 


l-cos(^+i?l)] 
"l-cosU-i^l)'] 


n-1 


y^  -   sin(n-u)  ^  (sinUii^2  '  sin  v^^)  +   (i&g  -  ^^)   sin  n  ^ 


u=l 


(28) 


where 


for   0  <  ;»  <  i»L 

2 
for  ^2  <  ^  <  It 


f  =  cos  n^  for  ^-j^  <  ^  <  ^ 


f  =  0 
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In  each  case,  by  meanB  of  a  few  terms  of  this  type,  the  fuac- 
tion  H  can  "be  expressed  and  the  pressures  can  therefore  he  com- 
puted as  a  stun.  The  computation  of  the  resulting  actions  involves 
easy  integrations.  In  order  to  simplify  this  part  of  the  computa- 
tion, Sohngen  gives  the  expression 


I  In  Ll-cos(i&-a)^  sin  n'^  d'^  =  -g"  cos  n  ^-  -  sin  na  + 

n-1     -  -. 

2  V  co3L(n-u)  i^  +vcl\    +  1  (cos  na  -  cos  nd)  ln[l-co8(^-a)] 
n  X  ^        n-  u         n 


The  use  of  these  expressions  will  he  clear  from  the  examples 
that  follow. 

14.  Examples  of  application.  -  Let  the  expressions  of  the  pre- 
ceding section  be  applied  to  the  determination  of  the  instantaneous 
pressures  corresponding  to  the  rotation  of  the  elevator;  that  is, 
it  is  assumed  that  the  forward  part  of  the  wing,  corresponding  to 
values  of  ^  between  0  and  ^q,   remains  immovable  while  the  rear 

part  rotates  rigidly  about  the  hinge  located  at  the  point  of  the 
abscissa  x  =  -  cos  i^q*  '^^  angle  of  the  elevator  is  denoted  by  3 
(positive  downward)  and  the  primes  denote  the  derivatives  3*  «  d3/ds, 
P"  =  d^p/ds^.  Hence, 


z  «  3(cos;^0  "  cos  ^ ) 
and  therefore,  from  equation  (l), 

w/v  =  P  +  3«  (cos  1^0  -  COS  ^) 
(p»+0"  cosi&o-P"  cos  ^)  sin^  d^  "Bo+B^  cos  ^  +  B2  cos  2^ 


^^ 


H  =  w+V 


^ 
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where 


(' 


1      P"        ^ 

Bq  =  V  ^  3  +  2(3'  cos  lit  Q  +  2  3"  +  -4-  cos  2i»q  J 


Bj_  =  -  V  (2P'  +  P"  cos  ^q) 


B. 


V  P" 
2  '-^- 


cos 


2d. 


For  each  of  the  three  terms  of  the  expression  of  H,  equation  (28) 
is  applied  with  n  =  0,  n  =  1,  and  n  =  2,  respectively,  and  with 
^1  =  ^0*  ^2  =  «•  MaMng  the  substitutions  yields 


Off 


linii^ 


H  d^ 


cos  iJ'  -  COS  1? 


1  (Bo  +  B-L  COS  ,>  +  B2  COS  2iJ)  In  1  -  cos  {d   +  ^p)  ^ 

2  1  -  cos  (  1?  -  1^0) 


(n  -  ^o)  (^1  sin  ^  +  Bg  sin  2^^)  -  2  B2  sin  -^q   sin  iJ 
On  the  other  hand. 


An  =  - 


Pit 


0 


2(jt  -  ^o)  ?Rt 

w  d^  = (3  +p'  cos  ^q)  +  £11-  sin  ^ 
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Tliere  is  thus  finally  olitained  from  equation  (26) 

rt  f/2  PT^  =<  P(«-iJo)  +  P'Tsin  ^q  +   (n-^Q^  °°^  '^O   '  '^°'^   2  "^ 
r2P'(n-^Q)  +  P"(n-^o)  °°^  ^  ■*"  ?  ^^^  ^o]  ^^^  '^"  ■!"  ("-'^o)  ^^^  2*  + 

,  r  -,  o"!   i-co8(i3+^o) 

ll^g  +  2P'(cos  ,^0  -  cos  ^)  +  2  P"(cos  ^0  "  cos  iJ )  J  In  i.go3(^.^^) 


As  a  second  example,  the  pressures  in  the  case  of  the  stationary 
gust  are  computed;  that  is,  the  wing  is  assumed  displaced  with  the 
velocity,  which  is  constant  in  magnitude  and  direction,  encountering 
air  layers  that  move  in  a  vertical  direction  perpendicular  to  V 
with  velocities  that,  at  each  point  of  the  fixed  space,  axe  main- 
tained constant  in  time.  The  graph  of  w  along  the  wing  trajectory 
(shape  of  the  gust)  is  assumed  given  and  P  is  a  general  point  of 
the  chord  that  is  indicated  by  the  positions  corresponding  to  t  =  to 
and  t  =  tg  +  dt  in  figxire  6.  For  the  point  considered,  the  value 


x: 


of   I   w  dx  is  represented  at  the  first  instant  hy  the  obliquely 

hatched  area  and  in  the  succeeding  instant  "by  the  same  area  increased 
hy  the  horizontally  hatched  strip  and  decreased  by  the  vertically 
hatched  part.  The  intervening  variation  in  the  interval  considered 
will  be  represented,  except  for  infinitesimals  of  the  higher  order, 
by  the  quantity  (w^-w)  ds.  Hence, 


H  =  w  +  -v-     wdx  =  W2 


The  quantity  H  is  therefore  constant  for  all  points  of  the 
chord.  The  same  conclusions  evidently  hold  for  the  point  P  if 
the  point  P  has  not  yet  entered  the  gust.   (The  only  difference 
with  respect  to  the  preceding  case  is  that  in  this  case  w  =  0.) 
Hence,  from  eqixation  (26)  (the  expression  for  p),  the  only  nonzero 
term  will  be  cot  ■^/2.      It  is  therefore  concluded  that  for  the  wing 
that  crosses  a  stationary  gust,  whatever  the  form  of  the  gust,  the 
pressures  are  distributed  proportionally  to  cot  i^  /2.  The  same  result 
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could  also  be  deduced  by  considering  equation  (18),  for  which,  in 
this  case,  the  first  member  should  be  zero  from  the  hypothesis  that 
w  does  not  vary  in  time.   The  result  given  by  Klissner  that  the 
aerodynataic  actions  on  the  wing  that  enters  a  stationary  gust  have 
a  resultant  passing  through  the  focal  point  is  thereby  obtained. 
This  result  holds  for  the  case  where  the  values  of  w  do  not  vary 
locally;  in  general,  in  agitated  air  the  velocities  vary  rapidly 
with  time. 

By  following  the  analysis  of  the  problem  of  the  stationary  gust, 
the  case  of  the  elementary  gust  (the  step  diagram  in  fig-  7)  is  first 
considered.  D'uring  the  time  in  which  the  front  of  the  gust  lies 
within  the  wing. 


r> 


Ao=- 


n 


0 


w  dv  =  Aw 

Jt 


n 


n 


w  cos  i»  d*  =  2  slni^' 


Aw 


0 


_   ^'  -  sin  ii*'  , 
w  = Aw 


where  x'  =  -cos  {>'    is  the  abscissa  of  the  gust  front.   When  the 
entire  wing  is  enveloped  by  the  velocity  Aw, 


Aq  ^  2  Aw     Ai  =  0 
The  instantaneous  pressures  are  therefore  represented  by 


jtp  =  2pV  Aw  d  '  cot  -  in  the  first  case, 

2  ' 


and  by 


p  -  2pV  Aw  cot  -  in  the  second. 
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For  the  computation  of  eqimtion  (23).  only  during  the  phase  of 
the  crossing  of  the  gust  front  is  dw  different  from  zero,  and  there 
it  has  the  value 


dw  =  ^-  (1-cos  i3>  )  d^'  =  ^Vpi^  dx' 
It  n  ^  1-x' 

If  8  denotes  the  distance  of  the  midpoint  of  the  wing  at  the 
actual  position  from  the  gust  front,  for  which  the  values  of  p  are 
measured,  at  this  instant 


-n   p/pv  Aw  cot  ^  =  2 


-1 


RCs-xM'^^dx' 


The  pressure  majr  therefore  he  computed  as  a  function  of  E,  in 

the  case  of  the  gust,  "by  means  of  a  simple  Integration.   The  resultant 
pressures  may  therefore  he  expressed  in  the  form 

p  =z  p  +  p  =  PV  Aw  El  cot  2  (29) 


where  Ej  is  a  function  of  the  distance  S2_  =  s  +  1  of  the  front 

of  the  gust  from  the  leading  edge  of  the  wing.  The  function  is 
evaluated  in  tahle  2j  the  values  are  ohtained  from  reference  13. 
When  S2_  is  negative,  then  evidently  E2_  =  0. 

From  the  solution  relative  to  the  elementary  case,  the  solu- 
tion for  a  gust  of  any  shape  can  he  obtained  hy  suhstituting  in  eq^ua- 

pw 

tion  (29),  in  place  of  E^  Aw,  the  quantity   I   E]_  dw  (taken  as 

the  a\m.  of  the  products  of  all  the  variations  that  the  values  of  w 
undergoes  for  the  corresponding  values  of  E^)- 

15.  Profile  in  harmonic  oscillatory  motion.  -  The  same  relations 
permit  solving  the  case  of  harmonic  motion.  Assume  V  constant. 
Using  the  complex  variable  notation  yields 

w  =  We 
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where  ^  is  the  frequency  and  W  is  a  function  of  ^     but  not 
of  t.  Equation  (26)  for  this  case  yields 


n  p /apVe  *  «  cot 


C" 


n* 


d^  -  sin  ^ 


W(i»  '  )  di^ ' 
cos  ^ '  -  cos  -i^ 


0)  sin  ^ 

Jo 


di>' 


r>Tr' 


COS  ^ '  -  cos  ^ 


¥(x)  dx 


(30) 


where 


-.12 


cos  ^' 


For  w* 


V  «  We    =  We   (W  constant) 


(For  the  position  s  -  o , 


_       _    -a;a   rr  iOt    -'^'O 
wCs-a)  «  w(s)e    "We    •  e 


(Hence,  from  equation  (23"), 


Py/2pWe    cot  I  -oiJ     B(0)« 

Jo 


e<^adO 


(31) 


The  same  problem  can  be  attacked  by  making  use  of  the  general 
relations  Initially  given.  By  procedures  that  are  developed  in 
numerous  publications  and  that,  in  part,  are  herein  presented,  rela- 
tions are  arrived  at  that  are  equivalent  to  equations  (30)  and  (31). 
In  this  meinner,  which  is  more  rapid  than  direct  computation,  it  is 
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found  that  the  quajitity  that  enters  the  second  member  of  equa- 
tion (31),  a  function  of  the  reduced  frequency  u)"  =  Ql/2Y,  is  iden- 
tified with  the  quantity  denoted  "by  K    in  references  2  and  9  and_is 
therein  expressed  by  means  of  the  Hankel  function  of  parameter  oj  : 


\  = 


Ho 


(2) 


Ho 


(27 


iH- 


T2y 


Eeference  14  gives  a  tahulation  of  this  function,  which  is  rather 
important  in  the  study  of  the  aerodynamic  phenomenon  for  the  oscil- 
lating wing  and  which  is  reproduced  in  table  3.  Kussner  uses  instead 
the  function  T  correlated  with  K     by  the  relation 


T  =  1 


2,\ 


whereas  in  the  paper  "by  Kassner  euid  Fingado  (reference  15),  the  func- 
tion ?  =  1  -  A.  is  used  with  argument  V  =  ui/2. 

American  publications  use  the  function  C  introduced  by  Theo- 
dorsen,  ■which  has  the  same  definition  as  P. 

By  means  of  equations  (30)  and  (31),  the  pressures  are  easily 
computed  for  any  type  of  oscillation  (for  example,  translatory,  rota- 
tional of  the  entire  wing,  or  rotational  of  the  flap).  The  coeffi- 
cients of  the  aerodynamic  actions  have  been  determined  in  various 
publications.  A  complete  tabulation  for  the  case  of  a  wing  with  a 
flap  hinged  at  the  forward  edge  is  found  in  reference  16,  the  com- 
putations for  which  were  developed  by  the  national  institute  for 
theoretical  applications  on  the  basis  of  the  formulas  of  Kussner. 
In  a  recent  publication  by  Kussner  (reference  17),  the  case  of  the 
profile  with  a  flap  and  with  a  tab  hinged  to  the  flap  is  treated. 
The  fact  that  the  hinges  of  the  two  movable  parts  can  be  retracted 
with  respect  to  the  corresponding  leading  edges  is  taken  into  account. 

16.  Analysis  of  pressures  on  airfoil  in  motion  in  nonpertiarbed 
air.  -  Equations  (27)  and  (28)  are  considered,  with  the  assumption 
that  the  values  of  w  are  due  only  to  the  motion  of  the  wing  and  can 
therefore  be  expressed  by  equation  (l).  It  is  first  assumed  that  V 
is  constant.  Then 


Px 


H  =  w  +  - 


dw    ,  TT   ^Z        ^Z        1 


at^ 


dx  + 


a 

5t 


'X 


ox 
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Inasmuch  as  adding  a  constant  to  the  value  of  H  does  not 
change  the  result  of  equation  (26), 


; 


OX 


and  therefore 


H 


3i 


5t 


Px 


a2z 

dt2 


dx 


(32) 


The  first  term  in  the  second  member  together  with  the  corres- 
ponding term  contained  in  Aq  gives  rise  to  the  pressures  that  axe 
denoted  by  Pq : 


1V0, 


Pit 


p^/2pv2  =  cot  I 


5J 


/Mt 


diS>' 


sin  iS> 


dz      di?' 
5x  cos  i>'  - 


cos  i> 


(33) 


0 


where  bz/6x  is  computed  at  the  point  of  integration  x'  «  -  cos  ^'• 
These  pressures  are  those  that  are  obtained  if  the  wing  in  the  actual 
configuration  is  under  the  conditions  of  steady  motion. ° 

The  second  term  in  the  second  member  of  equation  (32)  with  the 
corresponding  term  contained  in  Aq  gives  the  pressures  that  are 
denoted  by  p-i  : 


itp-j^/SdV  =  cot  2 


bz 


diS'  -  2  sin  t» 


j: 


dz      d^' 
3t  cos  1^'  - 


cos  ^ 


(34) 


These  actions,  which  are  proportional  to  the  vertical  velocities 
of  the  points  of  the  wing,  have  the  characteristic  of  damping  forces 
and  can,  in  part,  be  interpreted  by  kinematic  considerations.  Thus, 
if  the  wing  is  displaced  without  rotation  with  vertical  velocity  v. 


In  fact,  substituting  in  this  expression  the  value  of  u/V  given 
by  equetion  (37)  of  reference  11  for  oz/by.     and  developing  the  computa- 
tion yields  the  value  ot    y    =   p/pV  expressed  by  equation  (31)  pre- 
sented herein. 
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the  intuitive  result  that  the  wing  is  in  the  same  condition  as  if  it 
were  at  an  angle  of  attack  v/v  is  obtained  from  equation  (34). 
Also,  in  regard  to  the  effect  of  a  torsional  motion,  a  qualitative 
interpretation  of  -p-^     can  be  given*.  If  the  wing  is,  for  example, 

in  diving  rotation,  it  behaves  with  respect  to  the  fluid  as  though 
it  were  cuinred  upward.  It  is  seen  from  such  considerations  that  the 
focal  moment  that  arises  from  this  effect  is  tiumed  in  the  opposite 
sense  to  the  angular  velocity  (damping  action).  This  consideration 
of  the  dynamic  curvature  has  been  put  at  the  basis  of  the  numerous 
approximate  investigations  on  the  aerodynamic  coefficients  of  the 
vibrating  wing  (reference  9).  These  considerations  would  lead  in 
substance  to  the  computation  of  the  values  of  p-]_  with  the  same  rela- 
tion (equation  (33))  that  holds  for  Pq,  in  which  Sz/V  dt  is  substi- 
tuted for  dz/dx.  This  procedure  leads  to  results  that  are  quantita- 
tively in  error  because  it  is  necessary  to  halve  the  second  term  in 
the  expression  of  p-j_. 

The  third  term  of  equation  (32)  gives  rise  to  the  pressures  P2 : 

Ox' 


P2/2P  =  sin  1^ 


d;J' 


COB  ;>  -  cos  d' 

0 


§  ^  (55) 
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These  pressures,  which  result  independently  of  the  vertical 
accelerations  of  the  points  of  the  wing,  represent  an  inertia  effect 
of  the  mass  of  the  circulating  air.  As  evident  from  equation  (35), 
the  pressures  P2  do  not  depend  directly  on  the  local  values  of  the 

accelerations,  but  on  the  entire  distribution.  Thus,  in  the  case  of 
the  oscillation  of  a  flap,  although  the  forward  part  of  the  wing 
remains  fixed,  there  are  pressures  over  the  entire  chord.  There  can 
therefore, be  no  distribution  of  masses  that  are  apparently  capable 
of  reproducing  the  Inertia  effect  of  the  medium. 

The  result  of  equation  (35)  is  of  interest  for  the  rigid  motion 
of  the  airfoil.  In  the  case  of  translational  oscillation,  pressures 
are  obtained  that  are  distributed  proportionately  to  sin  ^  and  give 
rise  to  the  same  resultant  as  though  the  mass  of  the  cylinder  of  air 
circumscribed  about  the  wing  underwent  the  motion  of  the  wing.  It 
is  easy  in  such  a  case  to  compute  also  the  actions  on  a  part  of  the 
chord.  In  the  case  of  the  rotational  oscillation  of  the  profile 
about  the  mean  point,  the  pressiires  of  inertia  still  correspond  to  a 
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system  of  masses  distributed  according  to  sin  ^,   but  the  total  mass 
in  this  case  is  equal  to  one-half  that  of  the  preceding  case.  The 
rigid  wing  thus  undergoes  in  its  motion  an  inertia  action  that  can  he 
thought  of  as  reduced  to  two  masses,  each  of  a  value  equal  to  one- 
half  of  the  mass  of  the  circvmiscribed  cylinder  of  air  concentrated  at 
a  distance  L/V32  from  either  side  of  the  mean  point  (always  in 
regard  to  the  computation  of  the  resultant  actions). 

It  is  observed  that  the  pressures  P2  have  values  that  are 

independent  of  the  velocity  of  advance  and  must  therefore  be  sustained 
in  air  at  rest.  For  such  conditions,  the  results  wovild  not  rigour- 
ously  apply  because  the  assumption  of  the  smallness  of  the  perturba- 
tions, as  compared  with  V,  does  not  hold.  The  results  nevertheless 
agree  with  those  that,  for  ajiy  particular  case,  have  been  obtained 
without  the  preceding  essuiiiptions  and  moreover  apparently  reproduce 
sufficiently  well  the  actual  phenomenon  as  it  is  found  from  some 
measurements  by  Cicala,  in  which  the  periods  of  the  oscillation  of 
a  wing  model  in  rarefied  air  and  at  normal  pressure  were  compared. 
In  the  case  of  the  phenomenon  of  the  wing  vibrations,  the  inertia 
pressures  are  not  of  great  importance,  whether  because  the  additional 
masses  represent  a  small  part  of  the  mass  of  the  structure  or  because 
the  previously  described  actions,  which  exist  independently  of  the 
velocity  V,  can  be  directly  included  in  the  computation  if  measure- 
ments of  the  mass  of  the  structure  are  made  by  dynamic  procedures. 

Equation  (35),  which  permits  computing  the  inertia  pressures 
in  closed  form,  is  derived  in  reference  2,  and  the  expression  for  p^^ 

is  contained  in  the  same  reference.  The  expression  for  pi  presents 
a  certain  difference  when  compared  with  equation  (34)  in  that  it  con- 
tains a  term  that  in  the  preceding  scheme  is  added  to  the  transitory 
pressxire. 

If  the  velocity  V  is  not  constant  in  the  expression  for  H, 
there  is  derived  from  the  term  -    ^  dx,  in  addition  to  the 

quantity  expressed  by  equation  (32),  the  term  jx    -A  dx,  which 

u  tj   ox 

because  H  is  defined  except  for  a  constant  is  written  z   dv/dt. 
Hence,  to  the  preceding  computed  pressures  there  are  added  the  pres- 
sures P3,  given  by 


•='    dt         \co8  d-    -   cos  d' 
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Prom  this  relation,  there  is  obtained,  for  example,  the  following 
result:  If  the  rectangular  wing  is  displaced  with  velocity  V  not 
constant  and  with  single  of  attack  a,  which  is  constant,  there  is  a 
force  normal  to  the  flight  path  represented  by  the  distribution 

dV 
P3  =  2pa^  Bin^ 

which  corresponds  to  the  mass  of  the  cylinder  of  air  circumscribed 
about  the  wing  subjected  to  the  acceleration  dv/dt. 

It  is  emphasized  that  the  pressures  f  =  Po+Pi+P2'*"P3  ^^®  still 
added  to  the  transitory  pressures,  depending  on  the  values  that  the 
quantity  w  has  assumed  in  the  preceding  instants.  If  the  values 
of  w  are  distributed  linearly  over  the  chord,  w  represents  the 
value  of  w  at  the  neutral  rear  point-  The  instantaneous  pressures 
have  the  resultant  passing  through  the  focal  point. 

The  decomposition  of  the  total  pressures  into  instantaneous  ejid 
transitory  pressures  has  no  absolute  character  in  the  sense  that  a 
part  of  the  instantaneous  pressures  can  be  combined  with  the  transi- 
tory (not  vice  versa,  because  a  term  containing  the  history  of  the 
motion  is  clearly  distinguishable  from  the  terms  depending  on  the 
actual  values).  It  nevertheless  appears  that  the  definition  given 
herein  is  more  natural  because  in  the  limiting  case  of  steady  motion 
the  actions  resulting  from  the  group  comprised  of  the  instantaneous 
and  the  trgmsitory  pressures  become  zero. 

17.  Remarks  on  treatment  of  unsteady  motion  of  wing  in  two- 
dimensional  case.  -  The  first  investigator  to  study  the  aerodynamic 
problem  of  the  oscillating  wing  was  Birnbaum,  who  made  use  of  the 
concept  indicated  in  Section  6  of  the  splitting  of  the  circulation 
about  the  wing  into  bound  and  free  components.  Equation  (7)  is 
integrated  for  harmonic  motion.  By  use  of  the  complex  variable 
notation,  in  this  case, 

7  =7  (x)  e^^'^* 
£  =  £  (x)  ei^t 
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As  is  easy  to  verify,  equation  (7)  is  obtained  if 


(i)  =  -u)e 


-<^xj'^  7(x'  ) 


e'^^c'  cLx' 


(36) 


where 


OJ 


iQ/V 


Inasmuch  as  the  free  circulation  is  zero  in  correspondence  with 
the  leading  edge  of  the  profile  (if  no  other  sources  of  vortices 
occvir  upstream  of  the  wing),  in  order  that  x  *  -1  and   £  =  0,  the 
lower  limit  of  the  integral  in  equation  (36)  must  he  eqxial  to  -1. 
For  the  points  of  the  wake,  the  integration  is  evidently  limited  to 
the  chord,  heceuse  "Y  is  zero  outside  the  wing. 

On  the  "basis  of  these  results  equation  (8)  assumes  the  form 


2jrwo  = 


•i: 


(T-^£)  dx' 
x'  -X 


"r  dx' 
_1  x'  -X 


Ul 


e  dx' 


-1 


glxDc"     yf^jr")    ix'i 


(jO 


-CiJX' 


n 


dx' 


X'  -X 


e"^^"    7(x")  dx" 


-1 


where  x'   and  r"  are  variables  of  integration  and 


W  =  WqB 


.jQt 


Bimbaum  takes  into  account  the  condition  of  separation  by 
expressing  Y  by  means  of  a  combination  of  functions  that  satisfy 
this  condition.  The  solution  is  sought  in  the  form  of  a  series  expan- 
sion in  povrers  of  the  reduced  frequency  Ql/2V.  The  series  converge 
ratiier  slowly  so  that  the  results  of  Bimbaum  are  applicable  only  to 
rather  low  ■^'alues  of  the  reduced  frequency. 


Wagner  (reference  IP)  considered  the  problem  of  unsteady  motion 
of  the  wing  of  infinite  span.  In  the  first  part  of  reference  12,  the 
treatment  refers  to  the  case  In  which  w  is  consteuat  over  the  entire 
chord  but  variable  in  time,  a  case  of  fundamental  importance,  as  has 
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"been  stated  in  Section  10.  Making  use  of  the  condition  of  separation 
that  imposes  the  circulation  of  the  wing  capable  of  giving  rise  to 
a  finite  velocity  at  the  trailing  edge,  Wagner  arrived  at  an  integral 
eq^uation  that  defines  the  circulation  in  the  wake  of  the  wing.  With 
a  procedure  "based  on  the  moment  theorems,  Wagner  gave  an  expression 
for  the  computation  of  the  lift  and  of  the  moment  on  the  wing  on  the 
basis  of  the  circulation  in  the  wake.  In  particular,  the  computa- 
tions were  developed  for  the  take-off  motion  of  the  wing.  The  case 
of  rotation  of  the  wing  about  a  point  of  the  chord  was  also  considered. 

The  problem  of  the  unsteady  motion  of  a  wing  was  also  considered 
by  Glauert.  In  an.   initial  paper  (reference  18),  he  made  use  of  the 
hypothesis  that  the  circulation  remains  constEint,  a  hypothesis  that 
considerably  limits  the  importance  of  the  results.^  In  a  succeeding 
paper  (reference  20),  he  took  into  account  the  variation  of  the 
circiiLatioa.  In  contrast  to  Bimbaum,  Glauert  sought  to  obtain 
directly  the  total  circulation  of  the  wing,  which  is  divided  into  a 
part  that  would  correspond  to  the  case  in  which  the  vortex  wake  would 
be  absent,  and  into  a  circulation  induced  by  the  vortices  downstream 
of  the  wing.  It  is  simple  to  show  the  equivalence  of  the  relations 
asstimed  by  Glauert  for  computing  the  lift  P  and  the  moment  M  with 
respect  to  the  center  point  of  the  wing  with  the  relations  that  are 
obtained  on  the  basis  of  the  bound  circulation.  From  equation  (4), 

P/PV  =     (g-£)  dx  =  K  -    €  dx 

(37) 


i: 


m/py  =  I   (g-e)  X  dx 


Considering  that  for     x  =  -1,    e  =  0,  and  for     x  =  1,  according 
to  equation  (3),     V£  =  -  dK/dt,   integration  by  parts  yields 

J-l  Vdt      J.i       S 

J^_^    cxdx-  -2^p-^-2j   _^x2  3^dx 


Lamb  (reference  19)  also  treated  the  problem  with  the  same 
restriction. 


38 


NACA  TM  1277 


With  the  aid  of  these  relations,  equations  (37)  "become 


^/"-^^i*^ 


m/p  »  V 


gx  dx 


•\ 


(38) 


,    1  dK   1  S 
"^^  ■"  2  dt  -  2  5t 


gx^  di 


J 


considering  that,  according  to  equation  (5);, 


and  that  the  sign  of  the  differentiation  can  be  taken  outside  the 
integrals  because  the  limits  are  independent  of  the  time.  Equa- 
tions (38),  which  are  the  equations  used  "by  Glauert,  serve  for  the 
computation  of  the  resulting  aerodynamic  actions  "but  do  not  lend 
themselves  directly  to  the  determination  of  the  press\ire  distribution, 
a  computation  that  is  necessary  if  it  is  desired  to  know  the  actions 
on  part  of  the  wing  (flap).  The  equations  obtained  by  Wagner  are 
also  subject  to  the  same  limitation. 

With  the  values  of  g  expressed  with  the  aid  of  the  instantaneous 
characteristics  of  the  rigid  motion  of  the  wing  on  the  basis  of  the 
circulation  existing  in  the  wake,  Glauert  determined  the  lift  and  the 
moment  with  the  aid  of  eq_uations  (3S)  and  arrived  at  expressions 
agreeing  with  those  of  Wagner.  For  harmonic  motion,  in  which  the 
circulation  is  distributed  in  the  wake  according  to  the  sinusoidal 
law,  the  solution  of  the  problem  reduces  to  the  determination  of 
certain  integrals  that  Glauert  obtained  by  approximate  procedures 
that  limit  the  resiilts  to  values  not  much  higher  than  the  reduced 
frequency. 

On  the  basis  of  the  work  of  Bim"baum,  Kiisensr  (reference  8)  again 
took  up  the  problem  of  the  oscillating  wing,  assuming  for  the  bound 
circiilation  the  functions  that  are  used  for  steady  motion  (refer- 
ence 11,  p.  184).  The  corresponding  w,  expressed  as  a  function  of  x, 
consists  of  a  polynomial  and  of  a  trigonometric  function  multiplied 
by  a  factor  depending  on  ^,   which  Kiissner  computed  by  means  of  a 
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power  series.  After  rather  laborious  computations,  Kiissner  obtained 
the  pressure  distribution  on  the  wing  corresponding  to  the  motion  of 
translational  or  rotational  oscillation  of  the  wing  and  also,  approx- 
imately to  the  motion  of  the  flap,  for  which  he  derived  the  coeffi- 
cients of  the  aerodynamic  actions  for  the  field  of  variation  of  uj, 
which  is  of  interest  for  the  phenomenon  of  wing  vibration. 

Connected  with  the  investigations  of  Wagner  and  Glauert  is  the 
theory  developed  by  Theodorsen  (reference  4),  in  which  there  are 
separately  determined  by  use  of  the  methods  of  conformal  transfor- 
mation the  potential  fiinction  on  the  wing  for  simple  fundamental 
motions  of  the  wing  with  flap.  From  the  potential  function,  Theodorsen 
determined  the  pressures  with  the  aid  of  the  Bernoulli's  equation 
generalized  for  nonsteady  motion  (reference  11,  p.  38) : 

P  +  l!  +  ^  =  0  (39) 

P  .  2    5t  ^ 

(where  cp  is  the  potential  function). 

In  the  computation  of  the  difference  in  the  pressure  between  the 
lower  and  upper  surfaces  of  the  wing,  it  is  necessary  to  consider, 
according  to  epilation  (39),  the  q^uantity  derived  from  the  variation 
with  time  of  the  difference  in  potential  existing  between  the  two  sur- 
faces of  the  wing.  This  difference  is  meastced  from  the  circulation 
of  the  velocity,  which  is  determined  by  following  a  path  (fig.  3)  that 
Joins  the  two  points  situated  on  the  opposite  surfaces  and  passes, 
always  in  the  proximity  of  the  wing,  through  the  forward  edge  of  the 


x: 


wing,  thus  obtaining  the  quantity   j   g  dx;  it  is  thus  seen  that 

the  second  term  in  the  second  member  of  equation  (2)  represents  the 
coi*responding  term  in  dcp/dt  in  equation  (39). 

Theodorsen' s  treatment  of  the  problem  represents  a  marked  advance 
with  respect  to  the  preceding  work  because  it  determines  the  distri- 
bution of  the  pressures  on  the  wing  and  hence,  in  contrast  to  the  work 
of  Glauert  and  Wagner,  permits  the  computation  of  all  the  coefficients 
of  the  aerodynamic  actions  for  the  wing  with  a  flap.  Also,  because 
the  integrals  with  which  the  effect  of  the  vortical  system  of  the  wake 
is  computed  are  solved  by  means  of  Hankel  functions,  all  restrictions 
on  the  value  of  the  reduced  frequency  are  thus  eliminated  by  use  of 
the  existing  tabulations.. 
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Independently  of  Theodorsen,  and  almost  simultaneously.  Cicala 
arrived  at  the  solution  of  the  problem  of  oscillatory  motion  of  a 
deformable  wing  with  aoi'   law  "by  the  same  method  followed  by  Bimbaum 
and  Kiissner.  In  reference  9,  it  is  proven  that  a  class  of  functions 
exists,  depending  on  the  reduced  freq^uency,  that  represents  the  dis- 
tribution of  the  bound  vortices,  which  correspond  to  the  velocity  w 
distributed  over  the  wing  by  a  particularly  simple  law.  The  first 
of  these  functions,  corresponding  to  constant  w  on  the  chord,  con- 
tains the  Hankel  functions  of  the  reduced  frequency;  the  others  are 
essentially  represented  by  the  y^^     of  equation  (14)  and  give  rise  to 

the  values  of  w^  in  equation  (12).  In  this  manner,  relations  were 
obtained  by  means  of  which  the  coefficients  of  the  Fourier  series  of 
the  bound  circulation  on  the  wing  can  be  computed  in  closed  fonn  as 
a  function  of  the  coefficients  of  the  series  for  w.  The  coefficients 
of  the  aerodynamic  actions  for  the  wing  with  flap  were  thus  computed. 
In  a  succeeding  report  (reference  2),  it  was  also  shown  how  the  pres- 
sures depending  on  the  second  power  of  the  reduced  frequency  (inertia 
pressures)  and  those  proportional  to  the  first  power  (pressures  p-,  ) 

could  be  computed  in  closed  form  without  developing  them  into  Fourier 
series. 

At  the  same  time,  Kiissner  arrived  at  the  solution  of  reference  9 
by  a  procedure  described  in  reference  21,  some  results  of  which  were 
anticipated  in  reference  22.  The  general  case  of  nonsteady  motion 
was  also  treated  in  reference  21,  where  the  discontinuous  motion  was 
studied  on  the  basis  of  the  solution  for  the  harmonic  motion  with  the 
Mae   of  the  Fourier  integral.  The  case  of  a  stationary  gust  was  studied 
and  tests  were  conducted  (reference  23)  confirming  the  result  obtained 
that  the  focal  moment  remains  zero  during  passage  through  the  gust. 
The  solution  of  Kiissner  was  used  by  Dietze  (reference  24)  to  compute 
the  resultant  of  the  actions  on  the  flap  (in  the  preceding  papers 
only  the  hinge  moments  were  computed);  it  was  also  used  by  Krall 
(reference  16)  to  elaborate,  with  the  aid  of  the  National  Institute 
for  Applied  Computations,  the  tables  of  the  aerodynamic  coefficients 
for  the  oscillating  wing,  and  was  used  by  Dietze  again  (reference  25) 
for  the  computation  of  the  coefficients  for  the  wing  with  a  flap  and 
a  tab  hinged  to  the  flap. 

KJassner  and  Fingado  (reference  15)  also  succeeded  in  computing 
the  actions  corresponding  to  the  oscillatory  rigid  motion  of  the  wing, 
making  use  of  the  Wagner's  expressions  by  which,  in  the  case  of  har- 
monic motion,  the  evaluation  of  the  integrals  relative  to  the  effects 
of  the  vortices  in  the  wake  was  investigated  with  the  aid  of  Hankel 
functions  by  Borb^ly  (reference  26).  With  the  aid  of  the  Wagner's 
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expressions,  Ellenterger  computed  the  resultant  actions  on  the  wing  for 
a  flap  rotating  according  to  a  general  function  of  time  (reference  27 ) • 

In  a  summarizing  note,  Jaeckel  (reference  28)  established  a  coor- 
dination "between  the  procedures  of  Glauert,  of  Lamh,  and  of  Birnbaum- 
Kiissner  for  the  solution  of  unsteady  motion  of  wings  and  considered 
also  the  case  of  the  wing  with  a  variable  chord.  Jaeckel  also  pub- 
lished a  systematic  derivation  (reference  29)  of  the  results  that  were 
given  with  rather  synthetic  justification  by  Kiissner.  A  clear  review 
of  the  theories  on  the  wing  in  unsteady  motion  was  given  by  Lyon 
(reference  30). 

The  resiilts  of  the  preceding  studies,  with  some  further  develop- 
ment, are  treated  by  von  Karman  and  Sears  (reference  7);  a  derivation 
procedixre  is  developed  that  presents  in  an  intuitive  form  the  mathemati- 
cal fund.amentals  of  the  investigation.   The  computation  Is  restricted  to 
the  determination  of  the  lift  and  the  aerodynamic  moment,  which  are 
computed  with  the  aid  of  the  expressions,  respectively, 

p/p  =  ^z  ri  xi  (40') 

M/P  =|^z   r^x\  Uo") 

The  second  member  of  the  first  expression  represents  the  rate  of 
variation  of  the  moment  of  circulation  of  the  system  measured  with 
respect  to  a  fixed  point;  in  an  interval  of  time  dt,  this  moment 
should  vary  as  the  bound  vortices  are  displaced  by  the  amount  V  dt, 
while  the  position  of  the  free  vortices  and  therefore  the  moment  with 
respect  to  the  fixed  axis  have  not  changed.  The  total  variation  will 
therefore  be  given  by  the  product  of  the  total  botmd  circulation  and 
V  dt  and  therefore,  when  multiplied  by  the  density  and  divided  by  dt, 
will  give  the  lift.  By  ana,logous  reasoning,  the  second  of  equa- 
tions (40)  is  verified.  Thus,  the  square  of  the  distance  from  the 
fixed  point  varies  by  the  amount  SVX  dt  for  the  bound  vortices  while 
it  remains  constant  for  the  free  vortices.  Hence,  the  variations  of 
the  second  member  of  equation  (40" )  will  be  represented  by  the  moment 
of  the  bound  circulation  or  the  aerodynamic  moment  except  for  the 
factor  p . 

From  these  relations,  von  Ka'rman  and  Sears  derived  expressions  for 
the  lift  and  the  aerodynamic  moment,  which  present  a  better  generali- 
zation tlian  those  of  Glauert  Inasmuch  as  they  can  also  be  applied  to 
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the  deformable  profile;  "but  the  relatione  still  do  not  permit  compu- 
tations on  parts  of  the  wing.   The  authors  divided  the  totaJ.  circiila- 
tion  into  a  part  that  woxild  be  obtained  in  the  absence  of  the  vortex 
wake  (so-called  quasi- stationary  system)  and  into  the  induction  of  the 
vortices  of  the  vake.  The  case  of  the  general  motion  of  the  profile 
and  of  a  stationary  gust  is  also  considered. 

Garrick  (references  31  and  32)  "brought  out  the  relation  (Laplace 
transform)  that  exists  between  the  Wagner  function,  which  gives  the 
circulation  by  the  elementary  discontinuity,  and  the  function  that 
holds  for  the  h8.rmonic  motion  and  proposed  the  use  of  approximate 
expressions  to  represent  the  function  of  Wagner,  which  would  be  use- 
ful in  the  analytic  solution  of  various  problems  of  unsteady  motion. 

Posslo,  making  use  of  the  acceleration  potential,  considered  the 
problem  of  the  discontinuous  motion  of  a  wing  (reference  10)  and  the 
case  of  the  stationary  gust.  • 

In  a  recent  pu"blication  (reference  17),  Kiissner  and  Schwarz  indi- 
cate the  relations  ■^^dth  which  the  pressures  on  the  profile  can  be 
computed  without  making  use  of  ?ourier  series  but  using  integrations. 
These  equations  were  applied  to  the  determination  of  the  aerodynamic 
coefficients  for  the  wing  with  flap  ajid  hinged  tab  on  the  flap,  when 
the  case  is  considered  in  which  the  hinges  are  set  back  relative  to 
the  leading  edge  of  the  moving  parts.  The  same  relation  for  the  com- 
putation of  the  pressures  on  the  oscillating  wing  was  derived  in  a 
different  manner  in  a  report  by  Schwarz  (reference  5),  in  which  a 
clear  and  rigorous  derivation  of  the  known  solution  of  the  aerodynamic 
problem  of  the  oscillating  wing  is  developed. 

In  reference  13,  Kiissner  gives  the  general  solution  of  unsteady 
motion  in  the  two-dimensional  case  and  the  functions  of  Wagner  and 
those  relative  to  the  aerodynamic  actions  produced  by  a  gust  are  com- 
puted with  greater  precision. 

The  general  case  of  unsteady  motion  has  also  been  treated  in  a 
report  by  Sohngen  (reference  S).  The  solution  is  put  in  substantially 
the  form  indicated  in  Section  13. 
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PAET  II.  -  WBIG  OF  FINITE  ASPECT  RATIO 

The  relations  that  connect  the  velocities  induced  by  the  vor- 
ticity  of  the  inductor  system  of  a  wing  of  finite  aspect  ratio  axe 
derived  herein  and  a  form  for  these  expressions  is  sought  that  lends 
itself  to  a  future  refinement  of  the  investigation  that,  up  to  the 
present,  has  heen  conducted  with  approximate  procedures.  An  evalua- 
tion is  then  given  of  the  various  approximations  that  have  "been  used 
in  the  theory  of  the  finite  wing  in  unsteady  motion. 

The  consideration  is  restricted  to  essentially  rectilinear 
motion  of  the  lifting  surface;  that  is,  it  is  assumed  that  the  velo- 
city of  the  points  of  the  wing  give  small  deviations  relative  to  a 
mean  value  V,  which  maintains  its  direction  unchanged  with  respect 
to  the  motionless  fluid  at  infinity.  It  is  assumed  that  V  is  small 
compared  with  the  velocity  of  sound.  The  origin  of  the  orthogonal 
axes  X,  y,  and  z  is  located  at  the  point  at  which  the  Induced  velo- 
city is  measured  (point  of  induction).  The  x-aols  is  taken  parallel 
to  y  and  in  opposite  direction;  the  y-axis  is  normal  to  the  x-axls 
and  is  contained  in  the  fundamental  plane  that,  during  the  motion,  is 
at  a  very  small  distance  from  the  points  of  the  lifting  wing,  which 
Is  assumed  to  "be  of  infinitesimal  thlctaiess  and  curvature;  the  z-axls 
is  perpendicular  to  the  x-  and  y-axes  and  directed  downwards.  The 
vortices  having  an  axis  parallel  to  the  x-axis  (longitudinal  vortices) 
and  to  the  y-axis  (transverse  vortices)  are  considered  positive  if 
turned  in  the  sense  that  carries  the  positive  directions  of  x  and  y 
on  z.  In  the  middle  section  of  the  wing  is  located  the  origin  of  the 
|-  and  Tj-axes  parallel  to  x  and  y.  The  semispan  of  the  wing  is 
denoted  by  h,  so  that  for  the  points  -h  <_  tj  <_  h,   lo  ^^^     TO  denote 
the  coordinates  of  the  point  of  Induction,  hence 


^0 


y  =  T  -  no 

Le'fc  ^n(T)  "be  the  equation  of  the  leading  edge;  ^p(ri)  the  equa- 
tion of  the  trailing  edge;  ^p-^n  =  L  the  chord;  Q  the  pulsation; 
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(1)     th©  iniagiiiaxy  factor     iQ/Vj      ijj    the  reduced  frequency    C1l/27; 


s. 


■^  the  integral  taken  over  the  chord  from  the  leading  to  the  trailing 
/It 
edge.  The  other  symbols  are  defined  in  the  text  or  in  Part  I.. 

18«  Inductor  elements  in  tridimensional  case.  -  Let  P  (fig.  8) 
he  a  point  of  the  wing  that  moTes  relative  to  the  fluid  with  vector 
velocity  v  of  magnitude  v  and  having  any  direction;  the  linear  ele- 
ment dZ  through  P  normal  to  x  supports  the  hound  vortlcity  T; 
that  is,  supports  the  aerodynamic  action  pvFdl.   In  the  interval 
dt  that  precedes  the  actual  Instant,  the  point  P  starting  from  P' 
is  displaced  "by  the  segment  da  =  v  dt.  The  total  inductor  is  then 
changed  in  that  the  hound  vortlcity  element  has  come  to  occupy  the 
position  1-2  from  the  position  3-4,  leaving  behind  it  free  vortices 
and  simultaneously  creating  the  two  longitudinal  elements  2-3  and  1-4 
of  eq^ual  circulation  T  .  There  is  thus  added  to  the  preexisting 
inductor  system  a  closed  vortex  element  1-2-3-4;  on  the  side  3-4,  which 
at  the  instajit  t  there  exists  the  vortlcity  element  of  intensity 
-dP  liberated  from  the  hoiuad  vortex,  it  may  he  assumed  that  there 

dT 
simultaneously  exist  the  vortlcity  F -  —  da,  which  existed  at  the 

time  t-dt  and  the  element  of  intensity  T    with  oppositely  directed 
sign  constituting  the  fourth  side  of  the  circuit  1-2-3-4. 

With  the  aid  of  the  formxxla  of  Blot,  it  is  found  by  simple  com- 
putation that  the  closed  vortex  element  Induces,  at  a  point  0  at 
distance  r  from  P,  a  velocity  that,  except  for  infinitesimals  of 
the  higher  order,  nay  be  written  as^ 

d^  -  -  r  dZ  da/4rtr^ 

Whatever  happened  in  the  interval  dt  occurred  in  all  the  pre- 
ceding time  starting  from  the  Instant  at  which  a  force  on  the  lifting 
element  has  arisen;  that  is,  since  an  element  of  the  bound  vortlcity 
was  created.  In  the  problems  that  ordinarily  present  themselves,  it 
is  assumed  that  v  always  has  the  same  direction  parallel  to  the 
X-axis,  or  more  precisely,  the  deviations  with  respect  to  this  direc- 
tion are  considered  sufficiently  small.  Let  do  and  dl  therefore 

D  p 

The  symbol  d   is  used  to  indicate  a  quantity  that,  divided  by 
the  product  AZ  Ao,  has  a  finite  limit  when  these  elements  approach  zero. 
With  the  assumed  signs,  w  is  directed  upward  and  is  therefore 
ne^tive . 
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■fae  parallel  to  x    and  j,   respectively.  The  symbol  o    denotes  the 
distance  measured  from  the  point  P  of  coordinates  x     and  y, 
parallel  to  x,   assuming  that  the  circulation  I'  on  the  element  of 
the  lifting  surface  considered  is  expressed  as  a  function  of  a. 
Summing  the  effects  of  the  closed  vortex  elements  created  in  the 
preceding  time,  the  velocity  at  0  is  ohtained: 


dv  ^  -  =^ 


dy  jr  do 


4it 


Via 


0 


/\/(x+0)2  +  y2^ 


(41) 


The  velocity  induced  "by  the  hound -vortex  element  and  hy  the 
system  generated  by  it  for  r  ^  0  is  identical  to  that  which  is 
obtained  on  the  basis  of  the  concept  of  the  acceleration  poten- 
tial and  which  is  defined  as  the  field  of  the  pressure  dipoles 
(reference  33). 

The  particular  case  of  steady  motion  (  T   constant)  is  con- 
sidered. Inasmuch  as 


'^  Ob 


da 


0  V(x+a)^  +  y^' 


1  - 


X  +  Oi 


r\Jix-H3^)^  +   y2 


(42) 


a  quantity  that  is  denoted  by  f{x+Oi,j),   then  from  equation  (41), 

-  4«  dw  =  rf(x,y)  dy  (43) 

where  dw  denotes  the  velocity  induced  by  the  elementary  horseshoe 
vortex  of  frontal  side  dy  and  circulation  F .  The  velocity  (fig.  9b) 

induced  by  the  semi  vortex  2 »  is  -  (l  -  ^    "N   and  therefore, 

^^         VxS  +  y2  j 
summing  the  velocity  corresponding  to  the  other  semivortex,  which 
has  the  opposite  sign  to  the  first,  there  is  obtained,  except  for 
infinitesimals  of  higher  order,  the  quantity 


^•i[K^-^-)]^-^[-?^-^) 
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Adding  to  this  q.uantity  the  velocity  corresponding  to  the 
frontal  segment  yields  the  value  of  dw  given  in  equation  (43). 
The  equivalence,  which  in  this  case  holds  between  the  inductor 
system  of  the  pressure  dipole  and  the  elementary  horseshoe  vortex, 
can  also  be  seen  by  considering  that  the  sides  parallel  to  y  of 
the  continuous  vortex  circuits  (indicated  by  the  small  circles  in 
fig.  9(a))  are  canceled  if  p  has  the  same  value  for  all.  This 
inductor  element  is  called  a  bivortex  and  a  section  of  it  closed 
in  front  and  behind  is  called  a  segment  of  a  bivortex. 

Making  use  of  equation  (42)  and  integrating  by  parts  yields, 
from  equation  (41), 


4«  dw  =  r  f (x,y)  dy  +  dy 


\J 


g  f (x+o,y)  d        (44) 
0 


The  operation  that  leads  from  equation  (41)  to  equation  (44) 
transforms  the  inductor  system  constituted  by  the  closed  circuits 
into  a  system  of  bivortices :  a  generating  bivortex  originating 
on  the  wing,  denoted  as  a  bound  bivortex,  and  a  row  of  free  vortices. 
In  a  distance  do  of  the  wake  of  the  bound  bivortex,  free  bivortices 
originate  the  circulation  of  which  is  represented  by  the  variation 
that  the  intensity  of  the  bound  bivortex  has  undergone  in  passing 

through  the  distance  dO;  that  is,  2L  ^.a. 


do 


For  harmonic  motion. 


r  -  Te^^ 


where  V   is  constant  and  Q  is  the  frequency  of  the  motion.  In 
this  case,  V  is  always  assumed  constant.  The  circulation  of 
element  12  at  the  time  in  which  the  point  P  of  the  wing  was  set 
back  with  respect  to  the  actual  position  by  the  amount  j  is  repre- 
sented by 

r(o)  -  Te-^^o/y 

that  is,  with  respect  to  the  actual  circulation  T,  a  lagging  phase 
shift  represented  by  iTb/v  =  uo .  Making  use  of  this  resiilt  and 
denoting  by  ^V  dy  the  velocity  induced  by  the  bound  bivortex  and 
by  the  system  of  the  wake  in  the  case  of  harmonic  motion  yield 
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"Nod 


4jtw 


e-'^=^do 


0 


V(x-K3)2   ^  y2^ 


(45) 


By  setting 


X  =  Q  x/v  J  =C1  y/v"  0    »    Qo/v 


e<iuation  (45)  can  "be  vritten  In  the  form 


4nw  =  u?  $ 


(46) 


where  ^      is  the  function  of  two  variahles  x  and  y  defined  "by  the 
relation 


r\<xi 


$    (x,y) 


e-i<3  da 


[(i+0)2  +  ys]  ^/^ 


(46a) 


For     X  =  0,  this  function  has  "been  computed  "by  Miiller  (refer- 
ence  34,  p.   36).     For     x  ^  0, 


$  i^fj)  "   e 


ix 


«>(0,y) 


e    du 
(u2+5F2)3/2 


With  the  aid  of  this  relation,  which  is  easily  obtained  from 
equation  (46a),  the  computation  is  reduced  to  the  tabulated  function 
and  to  an  integration  "between  finite  limits  that  can  he  carried  out 
hy  graphical  or  numerical  methods.  The  computation  of  $  was  carried 
out  for  y  =  2.  When 


$"=$«+  i$" 

there  are  plotted  as  abscissas  and  ordinates  in  figure  10  the  real 
and  imaginary  parts  $«  and  $ "  of  the  function  and  the  curre  with 
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the  values  of  x  drawn.  This  curve  serves  for  the  computation  of 
the  velocity  induced  by  the  pressure  dipole  at  the  points  of  a  sec- 
tion at  distance  y  =  2v/^"^  from  the  inductor  elements.  The  point 
X  =  0  is  at  the  origin  of  the  dipole.  The  points  that  are  found 
upstream  correspond  to  the  positive  values  of  x.  For  this  segment, 
tlie  values  of  w  decrease  continuously  hut  all  have  slightly  dif- 
ferent phases.  Proceeding  toward  the  downstream  region,  the  points 
shifted  back  with  respect  to  the  front  of  the  bivortex  generator 
(x  <  0)  encounter  velocities  that  vary  little  in  magnitude  but  rela- 
tively more  in  phase.  It  is  evident  that  at  a  great  distance  from 
the  origin  the  velocities  vary  on  the  parallel  considered  by  a  sinu- 
soidal law  and  therefore  the  representative  point  of  figure  10  for 
negative  values  of  x  increasing  in  absolute  value  will  approach  a 
certain  limiting  circle  with  center  at  the  origin. 

19.  Correction  of  divergent  expressions.  -  The  expressions  of 
the  preceding  section  are  sufficiently  well  adopted  for  the  computa- 
tion of  the  velocities  corresponding  to  the  inductor  system  of  a  wing. 
Special  methods  are  required,  however,  for  the  points  of  infinity  that 
are  presented  by  the  functions  under  the  integrs-l  sign  present. 

The  case  of  steady  motion  is  first  considered.  A  bound  vortex 
filament  AB  (fig.  11(a))  of  circulation  F  varying  from  point  to 
point  is  considered  and  with  this  filament  is  associated  the  system 
of  longitudinal  vortices  that  is  shed  from  the  points  of  AB  according 
to  the  condition  of  Helmholtz.  The  velocity  that  these  elements 
induced  at  a  point  0,  which  is  outside  the  vortex  wake  of  the  fila- 
ment AB,  caji  be  computed  by  means  of  the  relation 

pA 
-  4jtw  =     rf(x,y)  dy  (47) 

Jb 

where  x  and  y  are  the  coordinates  of  a  point  of  the  filament  AB 
at  which  F  is  the  circulation.  This  expression  reduces  the  inductor 
elements  to  a  system  of  bivortices  that  originate  on  the  filament  AB 
(scheme  of  fig.  11(b)  for  A— ^0).  The  demonstration  of  this  eq illa- 
tion is  given  in  reference  35.  An  intuitive  Justification  is  obtained 
from  figure  11  (b).  The  segments  of  the  dotted  vortices  in  the  limit 
reproduce  the  effect  of  the  generator  filament  AB,  whereas  the  other 
longitudinal  elements  constitute  the  system  of  the  marginal  vortices; 
each  pair  of  contlgous  elements  function  as  a  single  inductor  ele- 
ment of  circulation  equal  to  the  difference  of  the  circulations  of 
the  two  elements  that  compose  it. 
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This  expression  is  applied  to  an  indefinite  vortex  of  unit  cir- 
culation at  a  distance  Xq  from  the  point  of  induction.  If  Xq  is 
positive,  f  alvays  remains  finite.  In  fact. 


2x' 


If  Xq  is  negative,  f  increases  indefinitely  for  y— K). 

Hence,  in  a  singular  integral,  it  is  necessary  in  computing  eq.ua- 
tion  (47)  to  exclude  from  the  integration  the  small  segment  from 
y=-5  to  y=+5  and  then  make  5  approach  zero.  The  expression 
that  is  thus  obtained,  however,  has  no  limit,  as  is  intiiltively  evi- 
dent from  an  examination  of  figure  12.  The  Mvortices  (fig.  11(a)), 
into  which  equation  (47)  transforms  the  inductor  system  when  the  ele- 
ments contained  in  strip  25  (which  includes  the  point  0)  are 
excluded,  are  equivalent  to  two  angular  vortices  (fig*  12(13))  (because 
the  semi vortices  Joined  "by  the  small  circles  are  canceled)  and  these 
elements,  on  approaching  0,  induce  a  velocity  that  increases  without 
limit.  For  this  reason,  the  integral  equation  (47)  diverges  if  the 
abscissa  xq     from  the  point  at  which  the  filament  AB  cuts  the  x-axls 

is  negative,  as  in  the  sinrplified  case  of  figure  12;  in  the  computa- 
tion of  the  principal  value  of  the  singular  integral,  an  expression 
is  obtained  that  increases  indefinitely  and  that  represents  the  velo- 
city induced  by  the  two  semi vortices,  which  are  indefinitely  removed 
from  the  point  0.   (See  fig.  12.) 

In  order  to  eliminate  the  previously  discussed  divergence,  the 
following  artifice  may  be  applied :  There  are  added  to  the  elements 
of  the  filament  AB  those  of  the  indefinite  vortex  I  (fig.  11) 
passing  through  the  point  Q  at  which  the  x-axis  through  the  point  0 
intersects  the  filament  AB  and  having  the  circulation  -Tq,  where 
Tq     is  the  value  of  P  at  Q.  The  system  thus  obtained  gives  a  velo- 
city at  0  that  is  denoted  by  Dw  and,  according  to  equation  (47), 
is  represented  by 


4n  Dw 


°[rof(xo,y)  -  rf(x,y)]dy  (48) 


The  first  term  of  this  integral  represents  the  effects  of  the 
bivortices  having  their  origin  on  vortex  I;  the  second  term  refers 
to  the  bivortices  having  origins  on  filament  AB  (hence  the  integra- 
tion is  extended  from  y  of  the  point  A  to  that  of  B).  The  two 
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terms  are  divergent  when  taken  separately,  "but  give  a  correct  expree- 
eion  when  added.  The  Intiiitive  reasoning  for  this  lies  in  the  fact 
that  the  two  semlvortices  that,  according  to  the  scheme  of  fig- 
ure 12,  include  the  point  0  and  axise  from  the  vortex  I  find 
compensating  terms  in  the  elements  arising  from  the  filament  AB. 
The  velocity  induced  by  the  vortex  I  is  then  added  to  Dw. 

In  the  case  of  a  system  Z  of  inductor  elements  analogous  to 
the  filament  AB  with  corresponding  marginal  vortices,  the  following 
procedure  may  "be  used : 

I)  The  transverse  vortices  cut  out  from  the  section  hy  the 
induced  points  are  prolonged  indefinitely,  thus  obtaining  a  system 
Ll>  the  induced  velocity  of  which  w*  is  computed  by  the  formulas 

of  the  two-dimensional  motion. 

II)  The  quantity  Dw,  which  constitutes  the  velocity  induced 
by  the  system  Z2  =  2  -  Z^  obtained  by  superposing  on  the  inductor 
system   £  the  system  Z^  vith  reversed  sign,  is  computed  with  the 
aid  of  equation  (47)  and  added  to  w' « 

The  advantage  of  the  preceding  procedure  lies  in  the  fact  that 
both  w'   and  Dw  have  a  homogeneous  inductor  system  (indefinite 
vortices  for  Zj,  bivortices  for  Zg)  ^^^   therefore  free  or  bound 

longitudinal  and  transverse  vortices  need  not  be  separately  considered. 
The  advantage  is  reflected  in  the  simplicity  of  the  formulation  of  w 
and  in  the  rapidity  with  which  the  results  of  the  existing  theories 
for  the  approximate  solution  of  the  problem  vmder  consideration  are 
derived  from  this  scheme. 

The  same  procedture  can  be  applied  if  the  motion  is  unsteady. 
The  bound-vortex  system  impresses  on  the  fluid  an  imprint  represented 
by  a  similar  vortex  configuration  corresponding  to  the  variation  of 
the  intensity  of  the  generating  system.  These  vortex  systems  carried 
by  the  stream  may  be  treated  as  shown  for  the  bound  system.  The 
plane  system  Z ,  including  the  free  vortices  can  be  analyzed  by  the 

formiilas  of  part  I.  The  system  Zg  is  made  up  of  bound  vortices 
with  the  corresponding  wake  of  the  row  of  free  elements  (that  is, 
pressure  dipoles)  and  can  therefore  be  analyzed  with  the  aid  of 
equation  (48). 

20.  Vortex  system  of  wing.  -  The  velocities  induced  by  the  vortex 
system  of  a  wing  can  thus  be  computed  by  adding  to  w' ,  induced  by 
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tlie  plane  system  that  would  ol^tain  if  the  traasverse  Tortices  cut 
out  from  the  section  through  the  point  of  induction  were  indefinitely- 
extended,  the  velocities  Dw  that,  as  a  function  of  the  "bound  vor- 
tices, are  written  as 

Dw  =    7w  di  dy  -     7q^  dx  dy  (49) 

Js        Jsi 

If  harmonic  motion  is  considered,  the  functions  w  can  "be 

derived  from  equations  (45)  and  (46).   If  the  intensity  of  the  gen- 
erator Taivortex  varies  according  to  a  different  law,  it  is  necessary 
to  make  use  of  equation  (41),  for  which  the  function  Tip)     must 
"be  known.  Prom  this  relation  the  function  ^  =  dw/dy  r(0)  is 
obtained,  which,  in  general,  varies  with  time. 

The  first  of  the  two  integrals  of  equation  (49)  is  extended  to 
the  surface  S  of  the  wing,  y     "being  the  "bound  vorticity,  the 
second  integral  is  extended  to  the  strip  S^  included  "between  the 

lines  Zi  and  t^     parallel  to  y  passing  through  the  points  at 

which  the  chord  through  the  point  0  cuts  the  leading  edge  I     and 
trailing  edge  t;  the  vorticity  -/q     on  this  strip  is  that  which  is 

on  the  wing  at  the  section  through  0.  On  each  element  of  the  sur- 
faces S  and  Sq_  there  originates  a  "bound  "bivortex  of  circulation 

7  dx  dy  connected  with  the  pressiire,  which  acts  on  the  element  "by 
the  relation 

p  =  PYy 

In  the  computation  of  Dw,  as  in  the  computation  of  w' ,  the 
total  circulation  g  along  y  can  "be  referred  to  instead  of  7, 
the  two  quantities  "being  connected  "by  the  relations  of  part  I.  When 
the  effects  of  the  elements  of  a  section  y  »  constant  are  considered, 
a  system  of  "bivortices  will  "be  obtained  the  origins  of  which  are  dis- 
tributed either  on  the  chord  or  in  the  wake,  the  intensity  of  the  dis- 
tri"bution  "being  represented  "by  g.  The  expressions  are  restricted  to 
harmonic  motion  and  are  written  as 


Q  A 

As  the  reduced  frequency  uu  increases,  the  modulus  of  w 
decreases  and  therefore  also  Dw,  which  represent  the  correction 
velocities  due  to  the  finite  span. 
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g  -  g(|^n)  e 


lOt 


■X 


Z  -    g  di  -  K(ti)  e 


iOt 


In  the  wake  at  a  distance   0  from  the  rear  edge  in  a  segment  do, 
■falTortlces  originate  that  have  the  circulation 


d  K(a) 
da 


do 


equal  to  the  change  that  the  total  circulation  K  has  undergone  in 
the  interval  in  which  the  rear  edge  has  passed  through  the  distance  dO 
at  the  time  t-o/v  in  which  Z(o)  was  displaced  with  respect  to  the 
actual  value  by  Qo  /v : 


K(o)  -  K  e-^^ 


The  velocity  induced  "by  the  elements  of  the  section  considered 
is  therefore  proportional  to 


Jgf  dx  + 
L 


OOO 


iOtel^f  do 


do 


gf  di 


Jo 


e-'^^f  do        (50) 


The  symbol     j_     denotes  the  ahscissa  of  the  rear  edge  of  the 
chord  (that  is,  Xp  ■  i-g  -ig)'   therefore 


f'(x,y)  -  f(x,y)   -  f(xp,y) 

It  is  easily  found  "by  using  equations  (44),  that  equation  (46) 
can  he  written  in  the  form 


(^2^  =,  f(x,y) 


-J: 


f  (x+o,y)  e-^^^  do 
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Hence,  equation  (50)  can  "be  transformed  to 


I" 


dx  -  0)  Kiji  (x-D,y) 


(50') 


The  eaqpreBsion  for  tlie  velocities     Dw    at  the  point    ^o'  ^0     ^^ 
therefore  obtained : 


4it  Dw  = 


J  -1)        J  L(ti  ) 


g(e,n)  f'(2,y)  dl- 


j: 


L(Tlo) 


g(bno)  f'U,y)  dt  - 


UJ 


•         $(xp,y)  • 


K:(tj)  dTi  +  uj2 


*(5Eo,y)  •  E(tiq)  dTi 


(51) 


where  Jw\    indicates  the  integral  taken  from  ^^^(ti)  and  ^-(ti), 
and  with  Xq  the  quantity   i  p(ilo)  "  ^0  ^A"* 

The  first  of  the  integrals  in  the  second  memher  of  equation  (51) 
represents  the  induction  of  a  system  of  segments  of  hivortices  having 
their  origin  on  the  points  of  the  wing  and  ending  on  the  rear  edge  of 
the  chord,  and  which  are  of  intensity  equal  to  the  total  vorticity  at 
the  point  of  origin.  The  second  integral  represents  the  effects  of 
the  analogous  inductor  system  existing  in  the  strip  S^.  The  third 

integral  expresses  the  induction  of  the  wake  system  constituted  hy  a 
system  of  pressure  dipoles  with  origin  on  the  trailing  edge  of  the 
wing.  The  fourth  integr-al  refers  to  the  elements  of  the  wake  having 
their  origin  on  the  line  t^^. 

If 

Hi 


G(e,Tl) 


■I 

en(^) 


gCesn)  dC 


(52) 
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the  first  two  terms  in  the  second  mem'ber  of  eqioation  (51),  which 
corresponds  to  the  velocity  denoted  ty  Wg,  can  te  transformed  by 

integration  by  parts  so  as  to  assume  the  form 


r>b 


4jrw2  = 


-b 


di 


G(e>n)  di 


L(Ti) 


dii 


G(e.TiQ)  dn 


L(t1o) 


where 


■tsjx^+j'^ 


It  is  then  observed  that;  in  regard  to  the  effects  of  the  ele- 
ments of  the  wake,  the  difficulties  arising  from  the  divergence  of 
equation  (47)  cannot  appear  because  all  the  points  of  the  wing  are 
located  upstream  of  the  line  from  which  the  system  of  the  wake 
originates.  Hence,  for  these  elements  the  decomposition  into  the 
systems  Zi  and  Ig  ca^  ^®  avoided  and  the  total  w  =  w'  +  Dw  can 
be  directly  computed.  The  velocity  denoted  by  W3  is  thus  obtained 
and  is  expressed  by 


4itw^  =  oo' 


-b 


If  wi  denotes  the  velocities  induced  by  the  plane  system  that 
is  obtained  by  indefinitely  prolonging  the  total  vorticities  cut  by 
the  section  containing  the  point  0;  that  is,  if 


SjtWn 


Jl(ti, 


eii  ,^0)  ^i 


) 


the  resultant  velocities  are  expressed  by  the  sum  w  »=  w,  +  w„  +  w„ , 


w, 


"2 
which 


In  general,  w^  must  represent  the  preponderant  part  of  w, 

facilitates  development  of  the  methods  of  iteration  for  the  computa- 
tion of  the  function  g  for  an  assigned  w.  The  velocities  W2, 

which  represent  the  part  whose  computation  presents  the  greatest 
analytical  complexities,  are  independent  of  the  reduced  frequency  and 
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are  therefore  determined  only  once  for  any  frequency  of  oscillation. 
The  values  of  W3,  which  depend  on  the  frequency  of  the  motion,  are 

expressed  by  simple  integrals  and  are,  moreover,  zero  for  all  the 
distributions  of  the  circulation  (equation  (11))  for  which  the 
integral  K  is  zero. 

21.  Approximate  theories.  -  A  first  approximate  treatment  of 
the  problem  of  unsteady  motion  of  the  finite  wing  was  developed  by 
Cicala  (reference  36).  The  principle  of  approximation  there  assumed 
finds  simple  formulation  if  reference  is  DKide  to  the  expressions  of 
the  preceding  section.  Although  use  is  made  of  the  exact  solution 
regarding  w' ,  in  the  computation  of  Dw  the  segments  of  the  blvor- 
tices  of  £2*  "to  which  the  first  two  terms  of  equation  (51)  correspond, 

are  neglected  and,  moreover,  it  is  assumed  that  on  every  chord  the 
velocity  Induced  by  the  system  of  the  wake  is  constant  eind  equal  to 
that  which  would  obtain  if  the  system  started  at  the  point  of  induc- 
tion. Hence,  referring  to  equation  (51),  in  addition  to  neglecting 
the  first  two  terms  of  the  second  member,  in  the  computation  of  the 
other  two  by   $,  a  function  of  the  abscissa  and  of  the  ordinate  of 
the  point  of  t  with  respect  to  0,  the  value  corresponding  to 
X  =  0  is  substituted.  As  a  consequence  of  the  first  approximation, 
Dw=0  if  K=0.  Hence,  if  g  is  represented  by  a  linear  combi- 
nation of  gn  of  equation  (11),  the  corresponding  velocities  are  the 
Wn  of  equation  (12)  and  the  pressures  are  obtained  from  the  7^^ 
defined  by  equation  (14).  Reciprocally,  if  the  vertical  velocities 
on  the  various  chords  can  be  represented  by  a  combination  of  Wq, 
the  corresponding  circulations  and  pressures  can  be  computed  on  the 
various  chords,  section  by  section,  as  though  the  motion  were  two 
dimensional.  This  fundamental  simplification  permits  reducing  the 
tridimensional  problem  to  a  single  case;  for  example,  that  of  ver- 
tical velocities  constant  on  the  different  chords:  any  distribution 
of  w  that  is  developed  in  the  series  of  equation  (16)  (in  this 
case,  Aji  is  a  function  of  the  coordinate  measured  along  the 
span)  requires  particular  examination  only  for  the  circulation 
corresponding  to  the  w  defined  by  equation  (17)  (this  simplifica- 
tion is  also, in  general,  variable  from  one  section  to  the  next). 

The  bound  circulation  corresponding  to  the  velocity 
w  =  w  =  constant  over  the  entire  chord  is  obt8.ined  on  the  basis  of 
eq-uaticn  (19)  in  which,  if  harmonic  motion  is  considered,  according 
to  equation  (22)  and  the  results  of  Section  15,  the  following  rela- 
tion must  be  substituted : 

C  =  w(l-A) 
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From  7,  C     is  ttien  conrputed  with  the  aid  of  equation  (36) 
and  then  the  total  circulation  yields 


K  '    I      (7+e)  dx 


X 


When  the  computations  are  made,  it  is  found  that  "between  w  _ 
and  K  or  between  the  total  amplitudes  of  the  c[uantitiee  W  and  K 
there  exists  the  relation 


O)   (H^^^O   -   iH(2)i)  ei'^  K  -  2WL 


where  H  is  the  Bankel  function  of  parameter  oa  . 

The  velocity  W  is  that  induced  "by  the  system  ^i;   to  this 
quantity  is  added  the  velocity  due  to  the  system  Sg,  \rtiich,  for  the 

assumed  simplifications,  is  also  constant  on  every  chord.  With  the 
computations  developed  in  reference  35,  equation  (51)  is  transformed 
into  the  relation 


i 


4jt  Dw  -   I        1^  N  dii  (53) 


where 


t]     coordinate  measured  from  middle  section  of  wing  parallel 
to  y  (which  is  measured  from  point  0) 

N  »  i  -  ooF 


F     function  of  variable     y  «  ^y/v     defined  by 

u-^y     <2^y2j 


.-iu/ 


du 


JO 
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This  function  is  computed  ty  this  relation  for  y  >  0.  For 
y  <  0,  F(y)  =  -F(-y).  The  function  F  is  tahulated  in  reference  37; 
in  a  recent  report  "by  Kiissner  (reference  33),  the  fuQction  N  is 
tabulated. 

Adding  the  velocities  induced  "by  the  systems  ^i  and  Zg  yields 

Oh 


^  ~  2L  "411 


-h 


H  N  dTi  (54) 


This  equation,  on  the  "basis  of  knovn  values  of  W,  defiiles  the 
distribution  of  the  total  circulations  along  the  spaji.  For  3=0, 
this  reduces  to  the  integrodifferential  equation  of  Prandtl.  Inasmuch 
as,  for  the  simplifications  assumed,  the  distrihution  of  y     over  the 
chords  is  similar  to  that  of  the  two-dimensional  motion,  on  the  "basis 
of  the  values  of  the  total  circulations,  the  problem  is  completely 
solved.  Given  the  series  of  equation  (16)  on  w,  the  circulations 
ajid  the  pressures  corresponding  to  the  part  that  is  developable  in 
the  series  of  vn  is  computed  as  if  the  motion  were  two-dimensional, 
whereas  for  the  remainder  w,  the  circulations  are  distributed  on  the 
chord  as  if  the  motion  were  two-dimensional  and  along  the  span  of  the 
basis  of  the  solution  of  equation  (54). 

In  a  succeeding  note  (reference  37),  the  procedure  was  applied 
to  the  determination  of  the  aerodynamic  coefficients  for  the  oscil- 
lating wing. 

Independently  of  reference  36,  Borb^ly  proposed  a  type  of  approxi- 
mation for  the  computation  of  w  for  the  finite  wing  (reference  38). 
Eeference  35  shows  that  the  expressions  that  Borb^ly  elaborated  for 
the  computation  of  Dw  in  the  particular  case  of  the  elliptical 
distribution  of  K  along  the  span  agree  with  the  results  that,  for 
this  case,  were  derived  on  the  basis  of  equation  (55). 

Possio,  concerned  with  the  problem  of  the  stability  of  small 
oscillations  of  the  wing  considered  as  a  rigid  body,  also  analyzed 
the  problem  of  the  oscillating  wing  of  finite  span.  Making  use  of 
the  concept  of  acceleration  potential,  he  derived  equations  (9)  and 
(10),  which,  however,  as  has  been  shown  in  Section  8,  are  also  justi- 
fiable on  the  basis  of  the  concepts  of  vortices.  The  solution  is 
expressed  by  Possio  in  the  form  of  a  power  series  of  the  parameter 
Qb/v".  The  value  limited  by  this  parameter  and  the  smallaess  of  the 
ratio  L/b  (large  aspect  ratios)  in  the  series  containing  the  powers 
of  ilL/2V  Justify,  for  the  computation  of  the  w*  corresponding  to 
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tlae  system  Zq^,  the  iatroductlon  of  slmplificationB  that  were  not 

adopted  In  references  36  and  37.  The  assumed  simplification  in  the 
computation  of  Dw,  when  it  is  reduced  to  the  scheme  represented  hy 
equation  (49),  can  he  thiie  defined:  The  function  w  is  computed 
by  means  of  equation  (46),  substituting  for  $   the  value  that, 
according  to  equation  (46a),  corresponds  to  x  -  0  and  that  is  there- 
fore constant  for  the  elements  of  each  chord.  Reference  39  contains 
the  principles  of  the  procedure.  Some  of  the  results  are  described 
in  reference  40  and  in  greater  detail  in  reference  41.  More  general 
cases  of  the  mbtion  are  considered  in  references  42  to  44.  Refer- 
ence 44  analyses  the  law  of  variation  of  the  lift  on  a  rigid  wing  of 
elliptical  plan  form  during  the  start  of  the  motion  (the  same  problem 
that  was  studied  by  Wagner  in  the  two-dimensional  field).  In  the 
computation  of  the  velocities  induced  by  the  transverse  vortices, 
there  was  assumed  (in  the  simplification  of  equation  (14)  of  refer- 
ence 44)  £in  approDcimation  different  from  that  used  in  reference  43; 
the  approximation  gives  for  the  case  considered  a  greater  precision 
of  the  results. 

Sears  (reference  45)  also  studied  the  problem  of  the  oscillating 
wing  of  rectangular  plan  form  with  approximate  procedures.  Reducing 
the  computations  to  the  scheme  of  the  preceding  section,  the  simpli- 
fication assumed  consists  of  the  suppression  of  W2,  while  a  rigorous 
computation  is  proposed  for  W3,  and  making  use  of  the  results  that 

are  obtained  for  the  infinite  wing  with  sinusoidal  distribution  of  the 
circulations  along  the  span.  The  computation,  which  is  intended  to 
eliminate  the  errors  inherent  in  the  approximate  theory  of  refer- 
ences 36  and  37  criticized  by  Sears,  does  not  give  results  more  accu- 
rate than  that  theory.  As  an  example ,_the  effects  of  the  inductor 
elements  of  a  chord  L,   for  the  case  U)  ■  1  ■'-*-',  are  considered  with 
the  point  of  induction  at  a  distance  y  ■  L.  Apart  from  factors  that 
need  not  be  considered  in  a  comparison,  the  induction  can  be  expressed 
by  the  equation  (50' ),  which  can  be  written  as 


y  Jz 


where  'I']^  is  the  nondimensional  factor 


In  the  cases  of  wing  vibrations  encountered  in  practice,  the 


reduced  frequency  turns  about  this  vsilue. 
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^l(x,7)=y2[f'(w)  -u,2$(^,"^j] 


^Ax,j)   =y^|  f'(x,y)  -0)^  $1-^,-^)1  -$'i  +  i$"i 


X  being  the  alasclssa  of  ttie  point,  with  the  vorticity  g  and  xp 
the  ahscissa  of  the  rear  end  of  the  chord  measured  vith  respect  to 
the  position  of  the  point  of  induction.  The  q.uantity  $i  vaB  com- 
puted for  the  elements  of  the  circulation  g  at  the  leading  edge, 
at  the  middle,  and  at  the  rear  edge  of  the  chord  considered,  and  in 
the  three  cases,  on  modification  of  the  position  of  the  induction 
point,  there  were  obtained  (see  fig.  13,  in  which  are  drawn  as  the 
ahscissas  and  ordinates  the  real  and  imaginary  parts  of  <i'i, 
respectively)  the  curves  I,  II,  and  III  in  the  figiire  with  the  values 
of  X  =  xp/L.  If  the  approximation  that  was  made  in  reference  36 
is  assumed,  in  the  three  cases,  for  any  position  of  the  induction 
point,  the  end  of  the  representative  vector  is  the  point  indicated 
"by  the  double  circle;  with  the  aid  of  the  approximation  of  Sears 
for  all  three  cases,  the  representative  points  are  those  of  curve  III. 
The  error  with  either  approximation  is  large.  If  the  greater  com- 
plexity of  the  computations  req^ulred  by  the  solution  of  Sears  is 
considered,  the  advantage  of  his  approximation  is  q[uestlonable. 

In  a  recent  publication  (reference  33),  Eiissner,  making  use 
of  the  acceleration  potential,  developed  a  new  approximate  theory 
of  the  oscillating  wing  of  finite  span.  The  approximation  assumed 
is  easily  related  to  the  expressions  of  the  preceding  section; 
referring  to  eq\iation  (49),  the  value  of  w  given  by  eq.ua tl on  (45), 
setting  u  =  X  +  a,  can  be  written  in  the  form 


4ffO»  e'^ 


\Jx 


e-^du 
V(u2  +  y2)3 


(45«) 


Kussner'  s  solution  can  be  obtained  by  setting  the  lower  limit 
of  this  integral  equal  to  zero.  The  Dw,  on  the  basis  of  this 

assimiptian,  are  proportional  to  e"^*  (where  i     is  the  coordinate 
measured  parallel  to  x  from  a  fixed  origin  of  the  wing  in  an  arbi- 
trary position);  that  is,  the  values  of  Dw  are  distributed  by  the 
sinusoidal  law  over  each  chord.  If  the  distribution  of  the  vorticity 
is  considered,  which  on  the  basis  of  the  solution  of  the  two-dimensional 
motion  corresponds  to  the  velocities  distributed  according  to  this 
law,  and  if  for  this  case  the  values  of  w  of  the  system  I>±    are 
summed,  velocities  are  obtained  that  can  be  represented  by  the 
expression 
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w  =  We      =  We  ^    ' 


(55) 


where  W  is  a  fuQction  of  r\. 

By  obtaining,  on  the  "basis  of  the  Bolution  of  the  two-dimensional 
prohlem,  the  velocities  induced  by  the  system  ^i     and  adding  the 
•values  of  Dw  computed  with  the  aid  of  eq^uation  (49),  simplified 
according  to  the  preceding  assumption  with  respect  to  the  limit  of 
the  integral  in  equation  (45)  and  transformed  hy  operations  analogous 
to  those  that  led  to  equation  (54)  (see  reference  35),  the  final  eq^ua- 
tion  is  obtained 


?? 


^(h(2). 


iH 


(2) 


2L(Jo-iJi) 


l)K' 


1_ 
4rt 


Pb 

dg' 

-b  "^^ 


N  dn 


(56) 


in  which  H  and  J  are  cylinder  functions  of  the  parameter  jo  a^^cL 


K« 


1 


ye 


'^Ul 


This  quantity  has  the  same  modulus  as  the  total  circulation  K, 
but  has  a  certain  phase  displacement  with  respect  to  it. 

By  solving  the  integrodifferential  equation  (56),  which  has  the 
same  kernel  as  equation  (54),  the  law  of  distribution  of  the  vorticity 
on  the  wing  corresponding  to  the  values  of  w  given  by  equation  (55) 
can  be  obtained.  If  it  is  assumed  that  i  -Vt  is  the  abscissa 
measured  in  a  system  of  reference  fixed  with  respect  to  the  fluid,  it 
is  concluded  that  equation  (55)  represents  a  distribution  of  velocity 
having  local  values  constajit  in  time,  as  would  be  the  case  of  dis- 
turbed air  that  presents,  along  the  trajectory  of  the  wing,  vertical 
currents  of  constant  velocities  in  time  (stationary  gust).  The  solu- 
tion indicated  refers  to  the  case  of  the  stationary  gust  of  sinusoidal 
form.  The  more  general  case  of  motion  can  be  studied  with  the  aid  of 
the  preceding  solution  and  the  solution  of  the  two-dimensional  problem 
when  it  is  considered  that,  according  to  the  approximation  of  Kussner 
(as  with  the  procedure  of  references  36  and  37),  the  values  of  Dw 
are  zero  for  the  distribution  of  the  g  that  gives  rise  to  zero 
values  of  the  integral  K    sjid  therefore  of  K' ,  which  has  the  same 
modulus  as  K.  Hence,  in  this  case,  to  the  values  of  w  represented 
by  a  combination  of  Wq  expressed  by  equation  (12)  there  correspond 
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the  values  of  g  expressed  by  equation  (11);  as  if  the  motion  were 
two-dimensional.  On  this  "basis,  reference  35  indicates  the  extension 
to  general  cases  of  the  solution  based  on  eq^uation  (56).  Eiissner 
makes  the  generalization  by  a  different  principle,  which  leads, 
however,  to  expressions  that,  in  the  limiting  case  of  the  infinite 
wing,  give  the  exact  solution  already  known. 

The  approximate  theory  of  Kiissner  therefore  cannot  be  derived  on  • 
the  basis  of  the  vortices  concept  developed  in  Section  18,  as  is  true 
for  the  solution  proposed  by  Cicala,  which  is  criticized  in  the  pre- 
ceding note  by  Eiissner  as  presenting  arbitrary  assumptions]  the  disa- 
greement between  equations  (54)  and  (56)  finds  its  justification  in 
the  different  principle  of  approximation  rather  than  in  a  fault  of  the 
derivation  method  based  on  the  vortex  concepts,  which,  according  to 
Kussner,  woiild  lead  to  erroneous  results.  The  approximations  thus  far 
assumed  all  lead  to  a  somewhat  inexact  value  of  the  induced  velocities, 
as  is  shown  in  reference  35;  the  theories  are  all,  except  that  of 
Sears,  constructed  so  as  to  converge  in  the  case  of  steady  motion  to 
the  theory  of  the  vortex  filament  of  Prandtl,  whose  approximation  has 
thus  far  been  proven  sufficient.  On  analytically  examining  some  local 
values  of  the  errors  conunitted,-'-!  all  are  shown  to  be  of  little  pre- 
cision, from  the  simplest  to  that  of  Sears,  which  consists  of  the  most 
laborious  application,  or  that  of  Eiissner,  which  is  based  on  the 
concept  of  pressure  dipoles.  Only  the  fact  that  in  a  limiting  case 

•'""'"If  the  solution  is  expressed  in  the  form  of.  a  decreasing  power 
series  of  the  aspect  ratio  A,,  it  is  found  that  the  error  in  the 
existing  theories  starts  from  the  term  in  logi\./A.  .  In  reference  35, 
rather  than  analyzing  the  order  of  magnitude  of  the  error,  it  is  pre- 
ferred to  carry  out  the  computation  for  concrete  cases  so  as  to  be 
able  to  compare  the  various  approximations. 

The  comparison  is  particularly  evident  by  making  reference  to  the 
concept  of  pressure  dipoles.  According  to  the  principle  followed  in 
reference  39,  the  dipoles  of  the  system  2  2  are  transported  parallel 
to  the  direction  of  the  x-axis  up  to  the  induction  pointj  according 
to  the  principle  adopted  in  reference  33,  these  dipoles  are  given  the 
same  displacment  ajid,  in  addition,  the  same  phase  displacment 
Q  x/v  (x  is  the  abscissa  of  the  dipole  with  respect  to  the  induc- 
tion point);  according  to  references  36  and  37,  these  elements  are 
given  the  same  displacement  and  phase  shift  Q(x-Xp)/v.  It  is  shown 
in  reference  35  that  the  three  approximations  alter  somewhat  the 
value  of  the  velocity  produced  by  the  pressure  dipole.  The  affinity 
of  the  three  principles  is  evident. 
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the  approximations  lead  to  a  theory  that  has  shown  itself  satis- 
factory in  applications  indicates  that  a  compensation  of  the  errors 
will,  in  a  certain  measure,  be  found  in  the  values  of  the  resultants 
of  the  actions o  A  solution  of  greater  rigor  would,  however,  he 
greatly  desirable. 
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PART  III.  -  DRAG  AM)  PE0PUL5IVE  FORCE^^ 

The  component  parallel  to  the  velocity  V  of  the  aerodynamic 
action  of  the  wing  of  infinite  aspect  ratio  in  nonsteady  motion  can 
he  readily  computed  on  the  hasis  of  the  solution  of  the  problem  of 
tvo-dlmensional  motion,  as  given  in  part  I.  Under  the  assumptions 
made,  for  any  law  of  motion  of  the  infinite  wing,  the  theory  per- 
mits computing  the  drag  or  propvLLsive  force.  For  the  hypothesis  of 
a  perfect  fluid,  the  profile  drag  is  not  considered,  nor  are  the 
variations  of  this  drag  due  to  the  unsteady  motion  computahle  with 
the  aid  of  this  analysis.  These  actions  are  therefore  added  to 
those  that  are  here  computed. 

For  the  wing  of  finite  aspect  ratio,  in  the  problem  under 
condide'ration,  the  uncertainties  mentioned  In  part  II  are  also 
encountered.  The  analysis  will  therefore  be  limited  to  the  results 
obtained  for  the  two-dimensional  motion. -^^ 

Drag  and  thrust  in  unsteady  two-dimensional  motion.  -  The 
symbol  E^  is  the  instantaneous  value  of  the  force  that  arises  in 
the  direction  of  V  on  a  segment  of  unit  chord  of  the  wing  in  a 
uniform  flow  of  velocity  V.  With  the  notation  of  the  preceding 
parts  and  in  the  same  range  of  validity  of  the  theory  there  given, 
Ri,  considered  positive  if  it  has  the  sign  of  a  resistance  and  neg- 
ative if  it  corresponds  to  a  propulsive  force,  can  be  computed  with 
the  aid  of  the  expression  given  by  Bimbaxim  (reference  3).-'-^ 

Ri  =  PV  I  ^  7  dx  -  JtPa2L/4  (57) 

where 

2a  =  lim^ ^^  7  sin  i^  (58) 


12 
The  numbers  of  the  figures,  the  equations,  and  the  paragraphs 

follow  froiD  the  preceding  part. 

•^^The  treatment  of  Schmeidler  (reference  46)  examines  the  aero- 
dynamic action  corresponding  to  assigned  vorticity  of  the  wing.  The 
method  cannot,  however,  be  generalized. 

The  sign  J^  indicates  the  integral  taken  over  the  chord  of  the 
wing  from  the  leading  to  the  trailing  edge. 
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The  aerodynamic  action  on  a  segment  dz  of  the  chord  is  repre- 
sented hy  the  force  PV7  dx  noznnal  to  the  line  of  the  axis  and  is 
therefore  inclined  to  the  z-axis  normal  to  V  hy  the  angle  dz/dx. 
The  integral  in  the  second  member  of  equation  (57)  therefore  repre- 
sents the  action  along  x  that  is  exerted  on  the  points  of  the 
chord.  The  negative  term  in  the  expression  corresponds  in  every 
case  to  a  propulsive  force  end  arises  from  the  suction  that  is 
exerted  at  the  leading  edge  hy  the  surroixnding  fluid  and  produces 
a  lowering  in  pressure  that  hecomes  infinite  for  the  wing  of  infini- 
tesimal thickness.  If  7  is  expressed  by  means  of  the  customsxy 
series  of  functions  cot  ^/S,  sin  ^,    ...,  sin  n^,   only  the  first 
term  can  give  rise  to  a  suction  at  the  leading  edge  because  the 
other  terms  represent  circulations  that  vanish.  The  quantity  a 
defined  by  equation  (58)  gives  the  coefficient  of  the  first  term  as 
found  ^Then  considering  that  for  1^ — >0, 

t» 

lim  sin  "^   cot  2  -  ^ 


The  steady  motion  that  is  considered  corresponds  to  a  chordwise 
distribution  of  vorticity  represented  by  the  same  7  as  for  the 
steady  motion.  Inasmuch  as  the  suction  at  the  leading  edge  depends 
on  the  instantaneous  value  of  7,  the  value  of  7  must  be  the 
same  for  nonsteady  motion  as  for  steady  motion.  In  addition,  when 
it  is  considered  that  for  steady  motion  the  resultant  force  in  the 
wing  direction  miist  be  zero. 


r 


0  =  PT 


(^)^7  dx  -  nPa^LA  (59) 


where  (dz/dx)o  is  the  slope  of  the  axis  of  the  wing  on  which,  under 
the  considerations  of  steady  motioa,  the  vorticity  distribution  7 
holds . 

If  Vy     denotes  the  velocity  induced  by  these  vortices,  then 


(s^j 


w^ 
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From  this  relation  and  from  eq_xiatlons   (57)  and  (59), 


Ei  =  P 


(v  ^  -  v^)  7  dx  (60) 


which  is  a  second  form  of  equation  (57)  given  hy  Jaeckel  (refer- 
ences 28  and  29)  and  is  entirely  eq.ulTalent  to  the  eq.uation  glyen 
in  the  development  of  the  computations. 

Denoting  "by  W2  the  velocity  induced  "by  the  free  vortices  and 
taking  account  of  equation  (1)  yield 

„  dz   dz 
Equation  (60),  on  the  "basis  of  this  relation,  hecomes 


Ej^  =  P    W27  dx 


i  7  dx  (61) 


This  third  form  of  equation  (57)  was  used  by  Schmeidler  (refer- 
ence 45).  It  lends  itself  to  interesting  interpretations:  If  the 
fox'ces  in  the  z -direction  are  distributed  with  density  PV7  and 
the  corresponding  velocities  of  the  points  of  application  are  3z/dt, 
the  instantaneous  power  N^     absorbed  by  the  motion  of  the  points  of 

the  profile  in  the  direction  normal  to  V  is  expressed  by  the 
relation 


WiA  =P  I  sir  dz  (62) 


It  is  noted  that  the  values  of  z  are  assumed  positive  down- 
ward, whereas  the  pressures  PV7  are  positive  upward.  Hence  W^, 

according  to  equation  (62),  is  positive  if  work  is  done  in  overcoming 
the  aerodyns/mic  action. 
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Wlien  equation  (62)  is  considered,  it  is  observed  that  the 
second  term  in  the  second  member  of  equation  (61)  represents  the  pro- 
pulsive force  (or  thrust )-'-^  that  vould  arise  if  the  phenomenon  occurred 
without  dissipation  of  energy;  that  is,  without  increase  in  the  kinetic 
energy  of  the  fluid  surroimding  the  wing.  The  first  term  represents 
the  dra^  R3  that  must  he  overcone  by  this  phenomenon,  or  in  other 
words,  by  the  creation  of  the  vortex  wake. 

The  computation  of  the  instantaneous  values  Hi  and  Ni  and, 
successively,  of  the  mean  values  Boi     and  Nja  of  the  same  magnltiide, 
in  the  csise  of  harmonic  motion,  is  immediately  obtained  on  the  basis 
of  the  expressions  given  in  part  I.  Poggi  (reference  47),  on  the 
basis  of  the  investigations  of  Glauert  (reference  20),  computed  the 
propulsion  and  power  corresponding  to  a  roteary  oscillation  and,  as  a 
limiting  case,  to  the  trans latory  oscillation  of  the  wing.  Kiissner 
(reference  21)  determined  these  values  for  the  wing  in  translational 
and  rotary  oscillation  of  the  wing  or  flap.  The  same  computations 
were  made  by  Garrick  (reference  48),  who  made  use  of  equation  (61) 
and  of  the  energy  interpretation  of  the  term 


Hm  =  |PLV2  1. 
2     \  nu) 


(63a) 


Hi(2)  +  iHo(2) 
-  Nm  =  IpLV^  *  (63b) 


where 


^^This  Is  the  quantity  called  "Vortrieb"  by  Schmeidler;  the 
first  term  is  called  "Widerstand . " 
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'2 


=  2i£.  f 

'    ^  Jl 


rtJj^  =  £iiil         Z  diJ  +  jtJ2 


rtJg  =  ^liH.  r   Z  cos  iJ   d^ 

and  z  =  Ze^'^  is  the  ordinate  of  a  point  of  the  wing  according  to 
the  complex  notation  with  ^(C),  |cL  and  C  the  real  part,  the 
modulus,  and  the  conjugate  complex,  respectively,  of  the  complex 
q.uantity  C,  with  H  the  Hankel  function  of  reduced  freq.uency 
0)  =Ql/2V. 

The  first  addend  that  appears  in  the  parentheses  of  the  expression 
of  I^  corresponds  to  the  resistance  Eg.  If  E  denotes  the  ampli- 
tude of  the  sinusoid  that  represents  the  distribution  of  the  free  Tor- 
ticity  in  the  wake,  it  is  found  that,  accoMing  to  equation  (63),  the 
following  relation  holds : 


Eg  =  pLE^/ie 


00 


This  relation  can  also  he  derived  "by  considering  that  the  work 
done  hy  the  resistance  Eg  during  the  displacement  2jt7/Q  corres- 
ponding to  an   entire  oscillation  of  the  wing  must  "be  equal  to  the 
increment  that  the  kinetic  energy  of  the  fluid  has  received  in  the 
same  time,  and  therefore  to  the  kinetic  energy  of  the  fluid  (considered 
stationary  at  infinity)  in  a  strip  included  between  two  parallels  to 
the  z-axis  at  2hY/C1   distance  from  each  other  and  located  in  the  wake 
at  a  great  distance  downstream  of  the  wing. 
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Applications  to  partlculsir  cases  of  oscillatory  motion.  In  the 
case  of  a  nondeforma'ble  wing  (a  ving  in  motion  of  trajislation  and 
rotation), 


Z  =  Zm  -  Zq  cos  1^ 


in  vhich  Z^  represents  the  complex  amplitude  of  the  oscillation  of 
the  middle  point  of  the  wing  ajid  Zq  the  amplitude  of  the  incident 
oscillation.  Equation  (63)  yields,  hy  simple  computations. 


vhere 


Ejj,  =  I  oLV2[5f(JZoK)    -  h   |J|2] 
-Njn  »  I  pLY^\jl  (  JZqK)   -  k  I  j|  2j 


(64a) 

(64b) 


J  =  Zo  +  IwZm  +  1.  iwZo 
2 


^  -        (2)''"'    (3)    <^-^-'>   ^  I  '^ 


Hi         +  iHo 


(2) 


^T^ni^ 


Hi(2) 


^°^L(2)        .      (2) 


Hi         +  iHo 
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The  quantities  K,  h,   and  k  are  functions  of  the  reduced 
frequency.  The  variation  of  h  and  k  as  f^^nctions  of  u;  are 
shown  in  figure  14.  In  figure  15,  the  real  part  is  plotted  on  the 
alDScissa  of  the  quantities  K/h  and  K/k  and  the  imaginary  part 
on  the  ordinate  assuming  the  segment  OA  as  unity  and  the  positive 
sign  of  the  imaginary  axis  in  the  downward  direction.  The  points 
on  the  curve  give  values  of  the  reduced  frequency. 

In  the  psirtlcular  case  of  translatory  oscillation  (Zo=0),  from 
equation  (64),  denoting  hy  v  »^|Zin|  the  maximum  velocity  corres- 
ponding to  this  oscillation, 

-%  =  I  PLv^ii  (65a) 

Nm  =  J  pLv^kV  (65h) 


The  force  that  arises  in  this  type  of  motion  is  a  propulsive  one. 
The  efficiency  of  the  ving  considered  as  a  means  of  propulsion  is 
h/k,  which  is  equal  to  1  for  ijj=0     and  decreases  continuously  toward 
0.5  as  u;  Increases. 

If  a  wing  of  velocity  V]^  normal  to  the  wing  velocity  V  is 

considered  under  conditions  of  steady  motion,  the  lift  P,  the  coeffi- 
cient of  which  is  equal  to  jtvi/V,  under  the  usual  assumptions  gives 
a  component  in  the  dif ection  V  represented  hy 

PV]^/Y  =  jtpLv^^ 


for  a  segment  of  unit  chord.  This  component  is  directed  forward 
whether  vi  is  directed  downward  or  upward.  On  varying  v,  harmo- 
nically, if  it  were  valid  to  apply  at  each  instant  the  expressions  of 
the  steady  motion,  vi  would  "be  the  mean  value  of  the  propulsive  force 
given  "by  equation  (65a),  in  which  h  =  1.  It  is  therefore  concluded 
that  the  exact  analysis  corrects  this  approximate  consideration  "by 
reducing  the  propulsive  force  hy  a  factor  depending  on  the  reduced 
frequency.  This  factor  approaches  1  when  the  reduced  frequency  is 
decreased  and  approaches  l/4  when  ou  is  decreased.  The  power  absorbed 
is  reduced  according  to  a  factor  that  varies  from  1  to  0.5  with  an 
increase  in  oj . 
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It  la  of  interest  to  examine  how,  by  combining  a  torsional  motion 
with  trans latlonal  and  oscillatory  motion,  the  propulsive  force  can  he 
Increased.  By  simple  computations,  it  is  found,  on  the  basis  of  equa- 
tion (64a),  that  for  a  given  amplitude  Zp  of  oscillation  of  the  rear 

neutral  point  of  the  wing,  for  every  reduced  oscillation  frequency  a 
certain  value  of  the  amplitude  and  of  the  phase  of  the  torsional  motion 
exist  for  vhlch  the  propulsive  force  is  a  maximum.  If  Z'  and  K" 
denote  the  real  and  imaginary  parts  of  the  quantity  K/h,  respectively, 
the  component  f  of  the  rotation  In  phase  with  the  trans latory  motion 
is  given  by 


zp  Ujg"  f... 


The  component  in  quadrature  is  expressed  by 


f£  2-K;.  00  (66b) 


The  propulsive  force  under  these  conditions  is  given  by 


^  =  f  p^'^-  •  ^4^  '«'' 


With  respect  to  the  propulsive  force  of  the  purely  trans latory 
oscillation  expressed  by  equation  (65a),  the  effect  of  the  rotation 
introduces  the  factor  dependent  on  the  reduced  frequency 
(K'2+e"2)/4(K' -l).  This  quantity  assumes  decreasing  values  with 
increases  in  the  reduced  frequency,  tending  asymptotically  to  the 
value  1,125.  For  U)   =  0.5,  this  value  is  equal  to  1.445;  for 
reduced  frequencies  not  too  small,  the  adding  of  the  torsional  motion 
does  not  greatly  modify  the  value  of  the  maximum  obtainable  propulsive 
force,  the  base  value  of  which  is  always  that  of  the  purely  flectional 
motion.  For  sufficiently  small  reduced  frequencies,  there  are  con- 
siderable Increases.  By  neglecting  the  higher  powers  of  the  reduced 
frequency,  the  expression  for  I^  can  be  assumed 


%  -  §■  Zp^V 
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cu,  the  torsional  motion  that  mast   he  comhined  with 


On  decreasing 

the  flectional  motion  to  ohtain  the  maximum  propulsive  force  tends  to 
assume  a  phase  displacement  of  90°  ahead  vith  respect  to  the  trans - 
latory  oscillation  (that  is,  in  the  phase  in  which  the  trans latory 
Telocity  is  a  maximum  upward,  the  wing  is  nearest  its  maximum  negative 
incidence  angle) . 


The  analysis  of  the  variation  of  the  propulsive  force  as  a  func- 
tion of  the  flectional  motion  for  a  given  amplitude  of  the  torsional 
motion  was  given  in  reference  49  by  means  of  a  graph  that  permits 
computing  directly  from  the  propulsive  force  (or  drag)  and  the  power 
ahsorbed  (or  emitted)  in  the  oscillation.  All  the  values  of  the 
efficiency  from  1  to  0  can  he  obtained  by  suitably  varying  the  ampli- 
tude and  the  phase  of  the  flectional  motion.  The  region  of  maximum 
efficiencies  is  nearest  the  point  J  =  0,  which  corresponds  to  the 
motion  without  drag  and  without  absorbed  power,  with  zero  vorticity 
in  the  wake.  Under  these  conditions,  for  a  small  ou,  the  torsional 
motion  is  displaced  by  about  a  90°  lag  with  respect  to  the  flec- 
tional motion;  that  is,  in  the  phase  in  which  the  trans latory  velo- 

city^^  upward  is  a  maximum,  the  wing  is  neejrest  the  maximum  positive 
incidence.  This  result,  in  relation  to  that  of  the  analysis  of  the 
maximum  propulsive  force  previously  indicated,  leads  to  the  conclu- 
sion that  the  conditions  of  maximum  efficiency  are  not  compatible 
with  those  of  maximum  efficiency  of  the  wing  considered  as  a  propul- 
sive means,  which  can  be  obtained  with  a  certain  loss  in  efficiency. 

On  the  basis  of  equation  (63),  it  may  also  be  determined  whether, 
by  combining  with  the  flectional  motion  a  deformation  that  alters  the 
curvature,  any  advantage  in  the  value  of  the  thrust  can  be  obtained. 


When 


then 


Z  =  Zm  +  Z2  cos  2^ 


J  =  2Z2  +  iooZj^ 


El 


m 


2^   JZ2 


-h  J 


•  I  PLT^ 


■^"^More  precisely,  the  velocity  of  the  rear  neutral  point. 
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In  this  case,  results  are  obtained  that  are  entirely  analogous 
to  those  of  the  motion  of  translation  and  rotation.  When 


2Z2  =  (f+iq)  Zm 


there  are  obtained  for  f  and  q.  the  same  expressions  of  equa- 
tion (66)  in  which  E'  and  K"  are  substituted  for  the  real  and 

imaginary  parts  of  Hi'^  Vh(Ei^^^+iHQ^^O  •  With  this  modification, 
the  factor  of  increase  in  the  maximvuu  thrust  has  the  same  expres- 
sion as  for  the  preceding  case.  This  factor,  which  can  he  expressed 
by  1/4 (k-h),  has  the  value  for  u) — ><*  and  the  value  1.148  for 
ou  =  0.5.  In  this  c£ise  also,  for  a  not  very  small  oj,  no  great 
increases  are  obtained  in  the  thrust  as  compared  with  the  purely 
translatory  motion.  As  in  the  preceding  case,  for  a  given  frequency 
of  oscillation  on  increasing  the  wing  velocity,  then  maximum  thrujst 
for  a  given  translatory  amplitude  first  increases  rather  slowly; 
only  when  sufficiently  low  values  of  the  ratio  Ql/2V  are  obtained 
does  the  maximum  thrust  tend  to  increase  linearly  with  the  velocity. 
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PART  rV.  -  EXPERIMENTAL  DTVESTIGATIONS 

The  experimental  investigations  that  have  thus  far  "been  conducted 
on  the  aerodynamic  actions  on  a  wing  in  unsteady  motion  are  not  as 
numerous  as  would  he  required  hy  the  complexity  and  Importance  of  the 
prohlem.  It  is  from  the  measurement  of  the  forces  on  the  oscillating 
wing  that  conclusive  data  are  expected  that  would  permit  a  reliable 
computation  of  the  critical  velocities  of  the  wings  and  tail  surfaces, 
Variovis  problems  relative  to  the  stresses  of  the  wing  structures 
during  flight  in  agitated  air  also  req^uire  experimental  clarification. 
The  experimental  investigation  should  furnish  the  necessary  control 
for  the  fundamental  hypotheses  of  the  theory  of  wings  of  infinite 
aspect  ratio  and  for  the  finite  wing,  the  actual  theory  that  makes 
use  of  approximations  that  have  not  yet  been  completely  checked  should 
be  integrated.  The  research  presents,  in  addition  to  the  difficulties 
common  to  all  problems  for  which  forces  variable  in  time  are  to  be 
measured,  serious  obstacles  for  the  requirement  of  absolute  regularity 
of  the  stream  in  which  the  experiment  is  conducted.  Small  fluctuations 
in  the  velocity  and  in  the  direction  of  the  wing,  which  do  not  have 
an^i'  great  effect  in  the  measurements  of  a  steady  flow,  can  render  the 
measurements  of  the  forces  on  the  oscillating  wing  entirely  unreliable. 

In  this  part,  the  results  obtained  up  to  the  present  by  various 
experimenters  will  be  discijssed,  and  the  results  compeired  with  theory. 

English  tests.  -  The  first  series  of  tests  was  conducted  by 
Duncan  at  the  National  Physical  Laboratory  and  published  in  1928 
(reference  50) .  The  object  of  the  tests  was  to  check  the  mechanical 
theory  of  the  wing  oscillations.  From  these  measurements  Duncan 
obtained,  for  a  particular  wing  model,  the  values  of  the  aerod5rnamic 
coefficients  to  be  introduced  in  the  expression  of  the  velocity  in 
order  to  compare  the  calculated  value  with  the  experimental  velocity 
obtained.  The  greater  part  of  the  tests  was  conducted  on  a  model 
that  was  deformed  during  the  oscillation,  according  to  an  incompletely 
defined  law.  The  tests  therefore  do  not  lend  themselves  to  a  check 
of  the  aerodynamic  theory,  a  check  with  which  the  experimenter  was 
not  concerned,  as  he  did  not  then  have  the  results  of  the  theory. 
A  series  of  tests  were,  however,  conducted  by  Duncan  on  a  model  that, 
during  the  oscillation,  rotated  rigidly  about  an  axis  parallel  to  the 
span.  The  wing  was  rectangular,  with  RAF  15  profile,  152 -millimeter 
chord,  and  686-millimeter  span.  The  axis  of  rotation  was  at  l/lO  of 
the  chord  from  the  leading  edge.  The  damping  of  the  oscillations  was 
measured  in  the  presence  of  a  wind  and  in  still  air  for  various  angles 
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of  attack  and  frequencies  of  oscillation.  Inasmuch  as  the  oscilla- 
tions did  not  have  very  rapid  damping,  the  results  can  "be  compared 
vlth  those  of  the  theory  on  harmonic  motion.  According  to  the  theory, 
the  moment  of  the  aierodynamic  force  due  to  the  rotational  oscillation 
possesses  a  component  in  phase  with  the  motion  and  a  component  ML 

in  quadrature  and  therefore  in  phase  with  the  angular  velocity  q; 
the  moment  Mg,  which  has  its  sign  opposite  to  q,  therefore  consti- 
tutes a  deunping  action  and  may  he  put  in  the  form 

Mq^  =  itphq  L^VS 

where  S  is  the  wing  area.  For  a  segment  of  an  infinite  wing,  the 
coefficient  b  depends  on  the  reduced  frequency  and  on  the  position 
of  the  axis  of  oscillation.  On  increasing  the  reduced  frequency,  h 
tends  to  the  valuel7 


40-1) 


(68) 


where  i     is  the  distance  of  the  axis  of  rotation  from  the  focus  con- 
sidered positive  if  the  axis  is  in  the  rear.  For  the  finite  wing,  h 
depends  also  on  the  plan  form.  Its  values  for  not-too-small  reduced 
frequencies  are  not,  however,  considerably  removed  from  that  given  by 
eqtiation  (68).  For  the  position  of  the  axis  of  rotation  of  the  tests 
of  Duncan,  this  gives  b  =  0.211.  The  values  of  b  obtained  on  the 
basis  of  the  damping  moments  measured  by  Duncan  are  given  in  figure  16. 
For  the  computation  of  b,  the  value  of  the  aerodynamic  damping 
moment  is  considered  to  be  the  difference  between  the  mesisured  values 
in  the  presence  of  wind  and  in  still  air.  The  values  are  all  below 
that  given  by  equation  (68).  For  equal  velocities,  these  values  indi- 
cate a  decrease  with  a  decreasing  u)^  as  would  also  be  given  by  theory 
for  positions  of  the  axis  of  rotation  ahead  of  the  focus.  For  equal 
frequency,  on  decreeising  the  velocity  (hence  on  increasing  jj) ,   the 
experimental  values  in  general,  indicate  a  decrease  that  can  be 
ascribed  to  the  effects  of  the  Reynolds  number.  There  is  also  a 
decrease  on  decreasing  the  angle  of  attack  (at  least  in  the  region 
investigated).  For  angles  of  attack  from  -4°  to  -5°,  Duncan  found 


■'-^Which  can  be  derived  on  the  basis  of  the  expressions  given  in 
references  2  and  21. 
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a  vanishing  of  the  aerodjnamic  damping.  This  phenomenon,  called  by 
Studer  "oscillations  of  separation,"  cannot  be  studied  "by  the  theory 
of  the  preceding  parts . 

The  measurements  conducted  by  Duncan  of  the  damping  of  the 
oscillations  of  the  flap  also  lend  themselves  to  a  comparison  with 
theory.  In  these  tests,  the  damping  due  to  the  friction  of  the 
suspension  is  rather  large.  By  assuming  for  this  case  also  that  the 
aerodynamic  damping  can  be  obtained  from  the  difference  between  that 
measured  with  wind  and  that  measured  in  still  air,  it  is  found  that 
the  experimental  value  of  this  damping  is  equal  to  about  one-'ialf  the 
theoretical.  This  disagreement  should  not  be  surprising,  because  the 
derivatives  relative  to  the  flap  are  always  markedly  less  than  the 
theoretical  values .  In  the  tests  conducted  by  Duncan,  the  value  of 
the  derivative  of  the  hinge  with  respect  to  the  angle  of  the  elevator 
under  conditions  of  steady  motion  was  eq^ual  to  0.6  of  the  theoretical 
value . 

Tests  conducted  at  the  Laboratorio  di  Aeronautioa  di  Torino.  - 
The  aerodynamic  actions  on  the  oscillating  wing  were  measured  by  Cicala 
in  the  free-Jet  wind  tunnel  of  600-millimeter  diameter  at  the 
Laboratoria  di  Aeronaut lea  di  Torino.  The  chord  of  the  models  on 
which  the  tests  were  conducted  was  about  13  centimeters  and  the  span 
about  50  centimeters .  Because  of  the  relatively  small  dimensions  of 
the  Jet,  which  was  free  in  the  region  in  which  the  model  was  located, 
the  wing  operated  with  a  rather  low  effective  aspect  ratio;  the 
value  of  5Cp/Sa  render  steady  conditions  (referred  to  pVS)  vas  eq.ual 

to  about  0.5rt  because  the  wing  projected  from  a  plane  that  was  placed 
tangent  to  the  Jet  in  order  to  mask  the  suspension  and  measuring 
apparatiis.  This  rather  low  value  of  the  aspect  ratio  is  one  dis- 
advantage of  these  tests,  which  are  described  in  references  14,  51, 
and  52. 

In  a  first  series  of  tests  (reference  14),  the  aerodynamic 
damping  of  the  flectional  oscillations  (rotation  of  the  model  about 
an  end  chord)  and  the  damping  of  the  oscillations  about  the  axis  con- 
taining the  foci  of  the  various  sections  (also  a  rigid  rotation)  were 
meas-ured.  The  measurements  relative  to  the  flectional  motion  were  of 
little  importance  because  of  imperfections  in  the  construction  of  the 
model  and  of  the  measuring  apparatus .  These  tests  were  later  repeated 
(reference  52).  The  measurements  relative  to  the  torsional  oscilla- 
tions gave  for  b  a  value  of  about  0.11  (against  0.125  given  by  eq.ua- 
tion  (68)  for  ^  =  0),  which  was  almost  constant  when  the  reduced 
frequency  was  varied  as  required  by  theory.  This  value  was  confirmed, 
at  least  for  the  range  of  not  very  large  angles  of  attack,  by  succeeding 
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tests  on  a  model  different  from  ttie  one  described  in  reference  14, 
this  model  also  being  of  symmetric  profile  but  of  greater  rigidity. 
The  principle  of  the  measuring  apparatus  for  these  tests  and  for 
those  of  series  III  can  be  briefly  described  as  follows:  The  oscil- 
lation of  the  model  was  controlled  through  an  intermediary  element 
that  possessed  two  simultajieouB  motions,  a  rotation  depending  on  the 
displacement  imposed  on  the  model  and  a  rotation  about  an  axis  per- 
pendicular to  the  first  rotation  and  depending  on  the  magnitude  of 
the  force  transmitted.  This  element  carried  a  mirror  that  reflected 
on  sensitive  paper  a  luminous  point  that,  by  describing  the  motion 
of  the  intermediary,  gave  the  force -displacement  diagram  (and  there- 
fore the  moment-rotation  diagram) .  The  test  conducted  for  equal  fre- 
quency in  still  air  and  in  the  presence  of  wind  permitted  isolating 
the  aerodynamic  action.  In  figure  17  are  given  some  of  the  oscillo- 
grams thus  obtained  that  give  the  simultaneous  values  of  the  angular 
position  of  the  wing  and  of  the  moment  transmitted.  The  field  of 
the  coordinates  can  be  retained  as  Cartesian,  so  that  the  diagrams 
are  approximately  ellipses.  The  enclosed  area  meaaures  the  work 
absorbed  in  the  oscillation.  In  the  figure  are  given  the  scale  of 
motions  and  also  the  lines  a  =  constant  conresponding  to  the  extreme 
positions  for  one  of  the  oscillograms  obtained  for  a  wing  velocity  of 
9.4  meters  per  second  and  for  a  number  of  oscillations  equal  to  570. 
For  each  measurement,  the  oscillogram  was  obtained  by  permitting  the 
luminous  point  to  rim  through  two  or  three  cycles .  The  paper  was 
successively  advanced,  thus  intercepting  the  light  point.  In  order 
to  obtain  a  reference  point  of  the  angles  of  attack  for  each  oscillo- 
gram, a  point  in  a  fixed  position  was  marked  (points  F  in  the  figure), 

In  the  third  series  of  tests  (reference  52),  the  damping  measure- 
ments of  the  trans latory  flectional  oscillation  (in  the  sense  pre- 
viously defined)  were  repeated.  The  component  of  the  lift  in  phase 
with  the  displacement  of  the  wing  can  be  expressed  by  the  derivative 


1   ^^P    „ 

where  v  is  the  trans lational  velocity  normal  to  V  and  v/v  is 
the  corresponding  variation  of  the  angle  of  attack.  According  to  the 
theory  of  the  infinite  wing,  the  factor  ag  varies  from  -1  to  -0.5 
on  increasing  the  reduced  frequency  (always  negative  because  the  veiri- 
ation  of  the  lift  has  a  sign  opposite  to  the  vertical  velocity). 
According  to  the  approximate  theory  of  the  finite  oscillating  wing, 
the  range  of  variation  of  this  factor  is  reduced,  starting  from 
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^  =  0,   from  the  correa ponding  value  of  the  steady  motion  computed 
vith  the  aid  of  the  vortex -filament  theory  of  Prandtl.  The  testa 
gave  for  ag  a  value  nearly  constant  and  equal  to  0.43  for  not 
too  small  a  velocity  within  the  range  of  reduced  frequency  in  which 
the  tests  were  conducted  (0.2  <  oj  <0.6).  At  small  velocities,  a 
decrease  of  ag  was  found  that  was  ascribed  to  the  effect  of  the 

low  Eeynolds  number. 

The  lift  component  in  phase  with  the  flectional  motion  and  the 
focal  moment  in  phase  with  the  torsional  motion  were  rather  small, 
as  required  "by  theory,  the  aerodynamic  inertia  effect  being  included 
with  the  measurements  of  the  mass  of  the  model  on  the  basis  of  the 
oscillation  data  in  still  air. 

The  variations  of  lift  that  arise  from  an  oscillation  of  the 
wing  about  the  focus  were  also  measured  (reference  51) .  By  expressing 
the  lift  component  in  phase  with  the  rotation  by  means  of  the  derivative 

1  SCp  ^  ^ 
jt  ha  ^ 


values  of  about  0.52  were  obtained  for  this  coefficient  for  all  the 
reduced  frequencies  at  which  the  tests  were  conducted  (0.1  <  i*J<0.7). 

The  component  Pq  of  the  lift,  proportional  to  the  angular  velo- 
city q  and  therefore  in  quadratiire  with  the  motion,  can  be  expressed 
by  means  of  the  relation 


Pg  =  rtpa4qSLY 


The  coefficient  a^,  which,  according  to  the  theory  of  the  wing 

of  infinite  aspect  ratio,  has  values  increasing  with  00  and  approaching 
to  0.5,  was  found  from  the  tests  to  be  almost  always  equal  to  0.38. 

The  principle  of  the  apparatus  for  the  measurement  of  the  lift 
due  to  the  torsional  moment  was  the  following:  The  wing  was  put  in 
forced  torsional  oscillation  by  guiding,  according  to  the  harmonic 
law,  a  point  of  an  end  section  while  a  point  of  the  saane  section  was 
attached  by  means  of  a  steel  wire.  Under  these  conditions,  in  addi- 
tion to  the  torsional  motion,  a  flexional  motion  of  rotation  arose 
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about  the  chord  of  this  section.  In  another  section,  a  force  having  a 
component  in  phase  and  one  in  q.uadrature  with  the  excited  motion,  was 
introduced,  the  amplitudes  of  which  could  be  varied  during  the  test. 
In  the  presence  of  wind,  the  amplitudes  of  the  two  components  were 
controlled  so  as  to  eliminate  the  flectional  motion  and  to  balance  the 
aerodynamic  action  and  the  inertia  force.  From  the  tests  conducted  at 
equal  frequency  in  still  air  and  with  wind,  the  aerodynamic  forces  were 
obtained.  By  the  same  principle,  the  focal  moment  due  to  the  flectional 
motion  was  measured.  In  agreement  with  theory,  the  value  of  the 
flectional  motion  was  so  small  that  It  could  not  be  measured. 

For  all  the  aerodynamic  derivatives  expressed  in  the  preceding 
form,  almost  constant  values  were  thus  obtained  in  these  tests  by 
varying  the  reduced  frequency.  The  theory  of  the  wing  of  finite  eispect 
ratio  developed  in  reference  37  justifies  this  result  for  the  range 
of  not-too-small  reduced  frequencies  and  low  aspect  ratios  at  -v^ich 
the  tests  were  conducted.  In  figure  16,  the  dotted  line  gives  the 
value  of  the  coefficient  b  that  would  be  obtained  by  these  tests, 
a  value  that  does  not  diverge  much  from  the  tests  of  Duncan. 

American  tests .  -  Tests  have  recently  been  conducted  in  the 
United  States  for  the  measurement  of  the  aerodynamic  forces  on  the 
oscillating  wing  to  check  the  theory  of  the  infinite  wing.  An 
interesting  series  of  testa  was  conducted  by  Held  and  Vincentl  at 
the  Guggenheim  Laboratory  (reference  53).  The  model  used  hsid  a  chord 
of  38  centimeters  and  therefore  permitted  the  attainment  of  suffi- 
ciently high  Reynolds  numbers.  The  span  was  not  large  (about  91  cm). 
Nevertheless,  a  large  aspect  ratio  was  obtained  because  the  model  was 
placed  between  two  walls  normal  to  the  plane  of  the  wing.  The  wing 
of  NACA  0015  profile  was  put  In  oscillation  about  an  axis  at  a  dis- 
tance of  4/10  chord  from  the  leading  edge.  At  the  opposite  edge  to 
that  at  which  the  motion  was  excited,  the  aerodynamic  action  was 
measured.  The  wing  support,  consisting  of  a  ball  bearing,  was  sus- 
tained by  a  rigid  spring  the  inflections  of  which  were  recorded  by 
means  of  mechanical  and  optical  amplification  on  a  strip  of  sensi- 
tive paper  with  uniform  forward  motion.  On  the  same  strip  were  marked 
the  instants  at  which  the  wing  occupied  the  extreme  and  middle  posi- 
tion. With  the  aid  of  a  harmonic  analysis  of  a  graph  of  the  forces, 
which  were  necessarily  irregular,  the  amplitude  and  the  phase  of  the 
fundamental  harmonic  with  respect  to  the  motion  of  the  wing  were 
derived  and  thiis  the  ratio  r  of  the  amplitude  of  the  lift  under 
conditions  of  oscillatory  motion  and  the  ratio  coirres ponding  to  the 
steady  motion  for  equal  rotation  and  phase  angle  (leading)   5 
between  the  lift  and  the  rotational  motion  were  measured.  The 
results  are  plotted  as  a  function  of  the  reduced  frequency  in 
figures  18  and  19  and  compared  with  the  theory  of  wings  of  infinite 


MCA  TM  1277  79 


aspect  rations  (continuoiis  curre)  and  with  those  that  vere  ohtained 
in  the  Torino  tests.  The  Torino  tests  must  necessarily  present  a  con- 
siderahle  divergence  from  the  American  tests  "because  of  the  difference 
in  aspect  ratio.  The  phase  displacements  predicted  by  theory  are 
somewhat  greater  than  the  experimental  values.  A  considerahle  devi- 
ation is  presented  by  the  theoretical  and  experimental  curves  of  the 
ratio  r. 

A  similar  series  of  testa  was  conducted  by  Silverstein  and 
Joyner  (reference  54).  The  model  hsid  a  chord  of  13  centimeters, 
considerably  smaller  than  that  vised  in  the  tests  by  Re  id.  In  this 
case  also,  a  large  aspect  ratio  was  attained  by  using  end  walls. 
The  oscillation  axis  passed  through  the  forward  quarter  chord  and 
the  aerodynamic  lift  force  was  measured  by  means  of  an  apparatus 
based  on  the  same  principle  as  previously  described.  In  these  tests, 
only  the  phase  displacement  between  the  lift  and  the  rotation  was 
measured  and  the  values  shown  in  figure  20  were  obtained;  the  con- 
tinuous lines  give  the  values  of  the  theory  of  the  infinite  wing  and 
the  dotted  line  give  the  values  obtained  from  the  Torino  tests.  The 
scatter  of  the  test  points  is  large  for  the  high  values  of  w ;  that 
is,  for  the  tests  conducted  at  low  velocity. 

Because  of  the  small  number  of  the  results  that  are  available, 
no  deductions  of  a  conclusive  character  can  be  given.  The  different 
conditions  under  which  the  tests  were  conducted  also  does  not  provide 
a  good  basis  for  comparing  the  different  results. 

The  tests  in  which  the  conditions  for  a  check  of  the  theory  of 
the  oscillating  wing  of  infinite  aspect  ratio  were  best  realized  are 
those  of  the  Guggenheim  laboratory.  The  comparison  is  not,  however, 
completely  satisfying.  The  probable  cause  of  the  divergence  encountered 
seems  to  lie  in  the  agglomeration  and  dissipation  of  the  wake  vortices, 
the  mutual  positions  and  intensities  of  which  the  theory  assumes  to 
be  maintained  indefinitely.  The  problem  should  be  investigated  more 

l^The  theoretical  curves  of  the  graphs  of  reference  53  do  not 
coincide  with  those  given  in  figures  18  and  19  because  the  aerodynamic 

inertia  effect  represented  by  the  terms  in  uu   in  the  expressions  of 
the  derivatives  is  not  considered.  In  fact,  this  action,  which  remains 
unchanged  with  and  without  wind,  is  already  compensated  in  the  pre- 
liminary operation  of  putting  the  center  of  gravity  on  the  axis  of 
rotation,  a  compensation  that,  it  seems,  was  effected  under  dynamic 
conditions .  The  correction  is  small  and  makes  the  test  points 
approach  the  theory  more  closely. 


80  NA-GA  TM  1277 


thoroughly,  especially  in  an  experimental  manner.  The  verification 
of  the  theory  of  the  infinite  wing  is  less  urgent,  however,  than  the 
investigation  of  the  finite  wing,  particularly  for  the  phenomenon  of 
wing  vibration,  in  which,  "because  the  motion  is  more  pronounced  toward 
the  tip  of  the  wing,  the  conditions  are  considerably  removed  from  two- 
dimensional  motion.  This  limiting  case  is  also  difficult  to  obtain 
experimentally  because  of  the  considerable  importance  assumed  by  the 
waJce  over  a  large  distance  behind  the  wing.  The  Torino  tests  make 
use,  however,  of  low  aspect  ratios  for  which  the  approximations  of 
the  theory  of  reference  31  are  less  Justified  for  giving  an  account 
of  the  results  of  such  tests .  There  would  therefore  be  required : 
First,  a  perfecting  of  the  theory  of  the  oscillating  wing  of  finite 
aspect  ratio;  and  second,  the  extension  of  tests  to  wings  of  greater 
aspect  ratio.  The  range  of  angle  of  attack  within  which  the  coeffi- 
cients can  be  held  constant  mijst  be  defined  and  the  field  of  coeffi- 
cients relative  to  the  oscillating  wing  with  flap  must  be  investigated. 


Translated  by  S.  Reiss 
National  Advisory  Committee 
for  Aeronautics. 
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